Quantum particles from classical statistics 
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Quantum particles and classical particles are described in a common setting of classical statistical 
physics. The property of a particle being "classical" or "quantum" ceases to be a basic conceptual 
difference. The dynamics differs, however, between quantum and classical particles. We describe 
position, motion and correlations of a quantum particle in terms of observables in a classical sta- 
tistical ensemble. On the other side, we also construct explicitly the quantum formalism with wave 
function and Hamiltonian for classical particles. For a suitable time evolution of the classical prob- 
abilities and a suitable choice of observables all features of a quantum particle in a potential can be 
derived from classical statistics, including interference and tunneling. Besides conceptual advances, 
the treatment of classical and quantum particles in a common formalism could lead to interesting 
cross-fertilization between classical statistics and quantum physics. 



I. INTRODUCTION 

Ever since the appearance of quantum physics it was 
widely believed that the basic concepts of a classical parti- 
cle and a quantum particle are fundamentally different and 
mutually exclusive. In this paper we argue that this is not 
the case. We describe a quantum particle in a setting of 
classical statistics with infinitely many degrees of freedom. 
On the other hand, we describe classical particles, with si- 
multaneously sharp position and momentum, within the 
formalism of quantum mechanics. Both quantum and clas- 
sical particles can therefore be described within the same 
conceptual setting. Their difference resides in the particu- 
lar dynamics, as expressed by different Hamiltonians which 
are based on different sets of observables. 

Perhaps the most striking evidence that classical par- 
ticles and quantum particles can be described within the 
same setting consists in the possibility of "zwitters" m - 
particles with properties interpolating between the quan- 
tum particle and the classical particle as a function of some 
continuous parameter 7. This parameter parameterizes dif- 
ferent possible laws for the time evolution of the probability 
density in phase space. While for 7 = a quantum particle 
shows interference in a double slit experiment, the classical 
particle for 7 = 7r/2 passes through only one of the slits 
without interference effects. Precision measurements can 
test the "quantumness" in a quantitative way by putting 
bounds on 7. 

Our implementation of a quantum particle within clas- 
sical statistics is not a deterministic local hidden variable 
theory but rather assumes that the fundamental descrip- 
tion of the real world is intrinsically probabilistic. We gen- 
eralize the description of quantum systems with a finite 
number of states M in terms of classical statistical ensem- 
bles [2-4J. The continuous density matrix and wave func- 
tion of a quantum particle is obtained by taking the limit 
M — 00. In the present paper we describe all observables 
of the quantum particle as standard classical observables 
that take fixed values for every state of a large classical 
statistical ensemble. Such an ensemble is described by a 
classical probability distribution. In two forthcoming pa- 



pers we will show that a much simpler description in terms 
of a probability density in usual phase space is possible if 
one admits the use of "statistical observables" . Statistical 
observables are computable in a non-linear way for a given 
probability density but cannot be associated to quantities 
taking a fixed value in every state. They are conceptually 
similar to entropy in thermodynamics. In this setting we 
can understand the quantum particle as a "coarse grain- 
ing" of a probability distribution of a classical particle. The 
present work provides the conceptual foundation for these 
further developments. 

The classical statistical foundation of quantum physics 
is based on several basic concepts 0-13]: 

(i) The quantum system is described as an isolated sub- 
system of a classical ensemble. This ensemble involves in- 
finitely many possible states for the particle and its envi- 
ronment. In a sense, the quantum particle is considered as 
an "excitation" of the "vacuum" (environment). 

(ii) The subsystem can be characterized by a restricted 
set of probabilistic observables, for which the spectrum of 
measurement values and the probabilities to find a given 
value can be computed from the information which specifies 
the state of the quantum system. 

(iii) In turn, the state of the quantum system can be 
characterized by the expectation values of a subset of these 
"quantum observables" . 

(iv) Conditional correlations, which are computable from 
the information contained in the quantum state the par- 
ticle, are used for predictions of the outcome of measure- 
ments of pairs of observables A and B. The conditional cor- 
relation or "measurement correlation" {AB)m differs from 
the classical or pointwise correlation {A - B). 

(v) In general, the joint probability of finding the value 
a for A and b for B cannot be computed from the quantum 
state. Even if joint probabilities exist, they typically in- 
volve information which characterizes the environment, in 
addition to the information characterizing the state of the 
quantum system. Concentrating on the quantum subsys- 
tem, we deal with "incomplete statistics" [H, where joint 
probabilities are not available or not used for {AB)m for all 
pairs of observables. We will find that a classical particle 
corresponds to the special case of complete statistics for po- 
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sition and momentum. Only for classical particles enough 
information is available such that joint probabilities can be 
used for {AB)m for all pairs of particle observables. 

(vi) Many different classical observables of the statistical 
ensemble, which describe the quantum system and its en- 
vironment, are mapped to the same probabilistic quantum 
observable. Their difference characterizes different proper- 
ties of the environment, while they give identical results for 
measurements concerning only the subsystem properties. 
The particle observables can be characterized by equiva- 
lence classes for the observables of the total classical sta- 
tistical ensemble. The measurement correlation {AB)m is 
a property of the equivalence classes, while the joint proba- 
bilities and the classical correlation {A ■ B) depend, in gen- 
eral, on the particular representatives of the equivalence 
classes or, in other words, on the environment. 

(vii) A time evolution which conserves certain statisti- 
cal quantities, as the purity and the copurity f^l of the 
subsystem, leads to the unitary time evolution which is 
characteristic for quantum mechanics. A more general time 
evolution can also describe decoherence Q or syncoherence 
i- 

In ref. [l]43| the basic conceptual settings have been de- 
scribed in detail for statistical ensembles which correspond 
to two-state and four-state quantum mechanics. It was 
shown how the quantum mechanical formalism with non- 
commuting operators arises form a description of a sub- 
system which obeys the basic concepts (i)-(vii). We have 
discussed explicitly quantum mechanically entangled states 
0) 0' B ^"^^ shown that the measurement correlation 
{AB)m is equivalent to the usual quantum correlation and 
violates Bell's inequalities 0. In this context the property 
of incomplete statistics 0, where the measurement corre- 
lation is not based on joint probabilities, is crucial 0,0- 
Complete statistical systems, for which the measurement 
correlation employs the joint probabilities, have to obey 
Bell's inequalities Usual local hidden variable theo- 
ries are assumed to be described by complete statistical 
systems. In contrast, for a fundamental probabilistic set- 
ting [l^ the assumption of a complete statistical system 
describing the subsystem can lead to severe problems 
This strongly suggests that incomplete statistics M is ap- 
propriate. Only in this case the EPR-paradoxon can be 
resolved satisfactorily and the entanglement, which 

is the basis for spectacular experiments on teleportation 
or quantum cryptography fl^l, or for quantum computing 
13| , can be explained. The use of probabilistic observables 
14| for the subsystem avoids conflicts with the Kochen- 
Specker theorem [l^, as demonstrated explicitly in 0. 

In the present paper we will not repeat the general dis- 
cussion of concepts and refer to [1| . Our purpose will be the 
generalization to quantum systems with a continuous fam- 
ily of quantum states, as appropriate for a quantum parti- 
cle whose wave function depends on a continuous variable 
as position or momentum. For this purpose we will first 
characterize the state of the particle-subsystem by a finite 
number of possibilities for the outcome of yes/no-questions. 
Such probabilistic systems, where the information concerns 
only Q bits, correspond to 2'5-state quantum systems. We 



will then take the limit Q — oo in order to construct the 
continuous location and momentum observables for a par- 
ticle. The start from a finite number of Q bits makes the 
possibility of incompleteness of the statistical description 
of the subsystem particularly apparent. 

In fact, the difference between a quantum particle and a 
classical particle precisely reflects the issue of incomplete- 
ness or completeness of the statistical description. For a 
quantum particle the joint probabilities for answers to loca- 
tion and motion questions are not available. This leads to 
non-commuting operators describing the position and the 
momentum of a particle and to Heisenberg's uncertainty re- 
lation. In contrast, a different set of yes/no-questions, for 
which joint probabilities for the outcome of location and 
motion questions can be given, describes a classical parti- 
cle. Now the position and momentum operators commute 
and a simultaneous sharp measurement for both types of 
observables becomes possible. For both the quantum par- 
ticle and the classical particle the formalism of quantum 
mechanics can be used. The description of the classical 
particle employs, however, an unusual set of observables 
where location and momentum operators commute. 

This paper is organized as follows: In sect. |ll]we discuss 
the position or location observable as an example of how 
to obtain observables with a continuous spectrum from the 
limit Q oo. In sect. Im] we then turn to alternative 
questions about motion and correlation for the particle. 
At this level the incompleteness characteristic for quan- 
tum particles plays a role. In sects. |IV]and|V]we discuss 
the observables describing the motion of the particle (mo- 
mentum) and the "particle-correlation", which is typical 
for a quantum mechanical wave packet. In sect. IVII we 
turn to the quantum mechanical description of a classical 
particle with commuting position and momentum opera- 
tors. We specify the Hamiltonian which leads to classical 
trajectories and show that the time evolution of the wave 
function implies the time evolution of the probability dis- 
tribution according to the Liouville equation for pointlike 
particles on classical trajectories. We compare the time 
evolution of quantum and classical particles in sect. IVIII 
Sect. I Villi deepens the connection between classical and 
quantum particles by using probabilities in classical phase 
space for the description of quantum particles. This is- 
sue will be discussed more extensively in two forthcoming 
papers. Conclusion and discussion follow in sect. IIXI 



II. PARTICLE LOCATION 

Consider a particle trapped in some region of space as, 
for example, a cavity or an atom trap. A priori we will 
not assume any property of the particle - it may be an ex- 
tended, pointlike of diffuse object, and one may be able to 
assign to it properties as momentum or angular momentum 
or not. It will be our task to find out under which condi- 
tions certain properties can be used for its characterization. 
We will assume, however, that the particle and its environ- 
ment can be described by a classical statistical ensemble. 
The environment includes the vacuum - which is a compli- 
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cated entity in the perspective of some more fundamental 
quantum field theory - and the trap. The classical statis- 
tical ensemble will be considered in the limit of infinitely 
many classical states and the classical ensemble is speci- 
fied by the probabilities for these states. One may view 
the particle as some type of excitation of the vacuum in 
presence of the "external fields " (typically electromagnetic 
and gravitational) which assure the confinement within the 
trap. 

In order to speak about a "particle" we consider our ob- 
ject as an "isolated subsystem" of the classical statistical 
ensemble. Isolation is used here in the sense that the ex- 
pectation values for a finite number of "basis observables" 
are sufficient for the characterization of the state of the 
particle. All properties of the particle should be expressed 
in terms of observables whose expectation values and cor- 
relations can be computed once the state of the particle is 
given. The state of the particle requires only a small part 
of the information contained in the probability distribution 
characterizing the ensemble of particle and environment - 
this is the basic reason why we typically have to deal with 
"incomplete statistics" Q, Q, 

Assume first that the properties of the particle reflect 
possible answers to a finite number Q of yes/no questions 
or Q bits. For example, we may use Q bits in order to 
characterize the location of the particle. For this purpose 
we ask Q questions of the type: "is the particle in a given 
portion of the volume?" , in order to divide the volume of 
the trap into M = 2^ cells. We enumerate these cells by 
a, a = 1 . . . M . The knowledge of the state of the subsys- 
tem requires then M probabilities > 0, X^q^" = 1' 
for the particle being found in the cell a. We may define 
M two-level observables A*^"^ that take the value 1 if the 
particle is in cell a and -1 if not, such that (A^"))^ = 1 
for all states of the ensemble. This induces for each cell a 
particle number A^'") = i(l + A^"^) which takes values 1 
or 0. 

We may further associate a cartesian coordinate Xi{a) to 
each cell - for example Xi{a) = /y^^-, UiCpy/ Jv(^a) '^i^h 
an integration over the volume V{a) of cell a. This allows 
us to define an observable for the location of the particle 
by 

A, = ^x,(a)A^("\ (1) 

with expectation value 

(A,) - X, = ^a:,(a)(Ar(«)) = (2) 

a a 

Here (A^^")) is computed in the classical statistical ensem- 
ble of particle and environment. Since A^^"^ equals the 
unit observable for a particle confined in the volume of the 
trap, one has indeed < (A^(")) < 1, = 1, as 

required for the identification {N'^"^) — Wa- 

We may also define observables for higher moments of 
the location observable, as the "dispersion tensor" 

Aff = ^ ia)x, (a) A^(") - x,x, , (3) 



with 

i^if) = (^'*^") ^ ^0 ~ ^j)Wa,. (4) 

a 

A "classical eigenstate" for the particle numbers consists 
of an ensemble for which (TV'")) = 1, {N^"^"^) = for 
some fixed index a. This implies also an eigenstate of Xi 

(2) 

with Xi = Xi{a) and vanishing dispersion (A]^ ) = 0. In 
contrast, for "equipartition" with {N'^°''>) = Wa — l/M and 
for cells with ^ii'^) — one finds Xi = and (A-\^ ') — 
For M — oo and fixed volume of the 
trap the location observable Xi becomes continuous. It is 
a possible candidate for a basis observable. 

We may cast our setting into a formalism familiar from 
quantum mechanics by associating the particle numbers 
A^*^"^ with operators 7V*^"\ which are defined as diagonal 
M X M matrices (no summation over a) 

(7V("))^^ = 6p^S^^. (5) 

We further introduce a hermitean density matrix p, whose 
diagonal elements read paa = Wa ■ One infers the quantum 
mechanical rule for the computation of expectation values 
in terms of the density matrix 

(7V(")) = tr(pA(")) =p„„ =7«„, 

(A,) = tr(pA,) , A, = ^x,(a)7V("). (6) 

a 

At this stage all operators are diagonal and only the diag- 
onal elements of the density matrix matter. There is no 
difference between "classical particles" and "quantum par- 
ticles" up to now. This difference arises only once we ask 
the question "What else can we know besides the location 
of the particle?" 

III. QUANTUM AND CLASSICAL PARTICLE 

In order to investigate the difference between quantum 
particles and classical particles we simplify our system to 
be one dimensional, choosing a trap in the form of a ring. 
The cell coordinates x{a) and the observable A are now 
periodic in the range — tt < x{a) < tt. For the location we 
can use a hierarchical sequence of bits, the first dividing the 
ring into its left and right half, the next for subdividing 
each half into quarters and so on. We will now consider 
more general types of questions. 

1. One-bit-particle 

Let us first investigate a hypothetical "particle" that can 
be described by only one bit, Q = 1, M = 2. By definition, 
it must be possible to characterize the state of this particle 
by only one yes/no question. This description of the sub- 
system by only one bit is considered as a basic property of 
the particle, while the probabilities w+ and W- for finding 
in an experiment the answers yes or no specify the state 
of the particle {w+ + W- = 1). The "particle observable" 
A is a probabilistic two level observable which can only 
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take the values ±1 (corresponding to the answers yes or 
no). In other words, its spectrum consists of two values 
7a = ±1. The probability w+ allows one to compute the 
expectation values for arbitrary powers of the observable. 
In particular (A^) = 1 is a sharp value independently of 
the state characterized by w^. The central issue is: what 
is the question? 

One possibility is the "location question", with associ- 
ated operators 

(7) 

However, other questions may also be asked as we demon- 
strate next in an example for a "one-bit-particle". (This 
"Gedankenexperiment" only illustrates the different possi- 
bilities a "one bit question" and should not be mistaken 
as a description of real measurements of a quantum parti- 
cle.) Suppose that we have placed in the ring two detec- 
tors on the right and left sides (say a,t x = ±7r/2). These 
detectors are assumed to signal if the "particle" goes 
through clockwise and —1 if it passes anticlockwise. We as- 
sume that the particle is not substantially disturbed by the 
countings of the detectors. Each one of the two detectors 
(l) (at X = — 7r/2) and (r) (at x = 7r/2) will record a series 
of -f 1 and —1 values. The probabilities to find given series 
of ±1 hits reflect the state of particle and environment, in- 
cluding the measurement apparatus. They contain much 
more information than contained in one yes/no question. 
From this information we want to extract the more limited 
information about the state of the particle using only one 
yes/no question. 

A measurement analyzes the series of ±1 values in a 
certain time interval At, which includes typically a large 
number of hits in the detectors. One possible two-level- 
observable for the particle takes the value +1 if there are 
more hits in (r) as compared to (l), and —1 in the oppo- 
site case. (For simplicity we assume an odd number of hits 
such that the question has a unique answer.) This "location 
question" can be associated with a "location observable" 
L. If there are more hits in the detector (/) than in (r) we 
define that the particle is in the left half of the ring with 
L = —1. The classical ensemble, which describes the par- 
ticle and its environment, specifies a probability for each 
possible series of hits. In turn, this determines the proba- 
bilities and wf^^ = 1 — to find for the location 
observable the values +1 or —1, which describe the state 
of the particle. If we associate to L the operator L = T3 
we recognize that X = {■jt/2)L corresponds to the location 
observable in different units. 

However, we may also ask a different question, as the 
"motion question": "do both detectors (r) and (1) to- 
gether show more -1-1 hits than —1 hits?". If yes, the "mo- 
tion observable" M takes the value (for clockwise mo- 
tion), if no, we associate M = — 1 to anticlockwise motion. 
Again, the classical ensemble induces probabilities ■u;^^'' 
and w'l*^'' = 1 — w^_:^^^ and therefore an expectation value 

for the motion observable (M) = w^'^ — wf^'' . The mo- 
tion observable can be associated with the "momentum" of 



the particle in appropriate units. The motion observable 
fulfills all criteria for a yes/no question for the one-bit- 
particle, just as the location observable does. 

A third two-level-observabic can be associated with the 
"correlation question" . We may register a "jump" if a hit 
in (/) is followed by a hit in (r), or if a hit in (l) occurs 
after a hit in (r). A "stay" is a sequence of two hits in the 
same detector. If a sequence of hits shows more jumps than 
stays we may assign the value -1-1 to the "correlation ob- 
servable" C, with the intuitive notion that for many jumps 
the particle has some tendency to be both on the left and 
the right side. In the opposite case of more stays than 
jumps the observable C assumes the value —1, with "anti- 
correlation" associated to the notion that the particle has 
a tendency to be exclusively either left or right. Again the 
average of the correlation observable can be expressed in 

terms of the probabilities w'^^ to find C = ±1, namely 

/r^\ (C) (C) 
(G) = w_^_ — w_ . 

The three two-level-observables L,M and C are unre- 
lated in the sense that the knowledge of the value of one of 
the observables does not yield any information on the other 
two observables. If we consider only chains with more hits 
in (r) than in (l), i.e. L = +1, this does not tell us anything 
on the relative number of -1-1 and —1 hits associated to M 
or on the relative number of jumps and stays associated 
to C. In other words, we assume that the restriction to 
ensembles with more hits in (r) than (l), i.e. an eigenstate 
to L with (L) = 1 , w^-* = 1 , wi^-* = 0, does not favor 
one of the two possible values for M or for C. Similarly the 
eigenstates of M and C should not "bias" the outcome of 
measurements of the two complementary variables. This 
can be realized for suitable rules for sequences of hits, for 
example by the "classical path sequences" where a jump 
consists of two hits of two hits —1 in the two different 
detectors (but not two hits with opposite sign) , while a stay 
consists of two hits with opposite sign in the same detec- 
tor (the particle moving forward and backwards through 
the detector). It also holds for "random sequences" where 
sequences with arbitrary and —1 hits are allowed. 

2. Particle subsystem 

On a somewhat more formal level we may associate each 
series of hits (each possible sequence of ±1 values) in the 
two detectors (r) and (Z) with a state t of the total sys- 
tem of particle and environment. The probabilities Pr 
for individual sequences characterize the classical statis- 
tical ensemble. They allow us to compute the probabilities 
w^\w'^^ and for finding the values -1-1 for the two 
level observables L, M and C, or equivalently to compute 
the expectation values (L), (M) and (C). We will now as- 
sume that the expectation values (L), (M), (C) character- 
ize the state of the one-bit-particle, and that this is the only 
information available for the subsystem. The joint prob- 
ability of finding L = 1 and M = 1 cannot be extracted 
from the knowledge of (L), (M), (C) - it is not a property of 
our one-bit particle. The joint probability would be com- 
putable from the knowledge of Pr for all sequences r, but 
this involves properties of the environment. 

In fact, if the joint probabilities for the dififerent combina- 
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tions (+, +), (+, — ), (— , +) and (— , — ) for the observables 
L and M would be available from the state of the par- 
ticle, we could construct a further observable LM which 
takes the value +1 if L and M have equal sign, and — 1 
for opposite signs. Its expectation value {LM) would be 
computable from the particle state. However, {LM) can- 
not be computed form (L) and (M) alone. Furthermore, 
the observable C is different from LM . We conclude that 
{LM) cannot be computed from the knowledge of {L) , (M) 
and (C), but would require additional "environmental in- 
formation" not available if the particle state is only charac- 
terized by (L), (M), (C). The statistical system describing 
the particle alone is incomplete. Furthermore, if LM would 
be a particle observable, then L, M and LM would be suf- 
ficient in order to construct a composite observable with 
a spectrum of four different values, corresponding to the 
(+,+),(+,—),(—,+) and (— ) combinations for L and 
M . A spectrum with four different values for a particle 
observable is in contradiction with the assumption that all 
particle observables use only one bit of information. 

Let us consider L, M and C as basis observables for the 
subsystem. The state of the subsystem is characterized 
by their expectation values (L) , (M) and (C). We may 
associate to these observables the operators (2x2 Pauli 
matrices) 

L = T3 , Af = -T2 , (7 = Tl. (8) 

We can further construct a hermitean matrix p which char- 
acterizes the state of the particle 

p=i(l + (L)T3-(M)T2 + (C)ri). (9) 

It obeys tip = 1 and < paa < 1- For all observables 
A = {L, M, C) the expectation values obey 

{A) = trip A). (10) 

The third condition for p being a density matrix, trp^ < 
1, is obeyed only for (L)^ + {M)"^ + {C)^ < 1 . This implies 
that at most one of the three observables can have a sharp 
value. For example, {M) = 1 implies {L) = (C) = 0. 
This is a typical situation for a quantum particle. For our 
one bit particle the sum (L)^ + (M)^ -I- (C)^ describes the 
purity P of the subsystem. We will assume P < 1, in 
accordance with the restriction to one-bit observables. In 
contrast, a classical particle would have a sharp value for 
both L and M and therefore (i)^ + {M)^ + (C)^ > 2. 
However, simultaneously sharp values of L and M answer 
simultaneously two yes/no questions ("is L positive?" and 
"is M positive?") such that this case should be described 
by a particle characterized by two bits. 

3. Two-bit-particle 

Let us therefore extend the discussion to a hypotheti- 
cal particle characterized by two bits, (5 = 2, M — A. 
We generalize the observables L,M,C to fifteen two- level 
observables A^'^) , fc = 1 . . . 15, {A'^^^)^ = 1. They are con- 
sidered to be the possible basis observables of the isolated 



sub-system, such that the state of the particle is described 
by fifteen real numbers pk , — 1 ^ Pfc ^ 1, 

Pk = (AW). (11) 
We define the purity of the subsystem as 

k 

and choose P < 3, as appropriate for two bits. Indeed, with 
Q yes/no questions we can fix at most 2*^ — 1 sharp values 
of independent two level observables, similar to the particle 
numbers A^(") in section ini (The constraint Y.a^''"^ = 1 
reduces the number of independent observables by one.) 

We introduce the M'^ — 1 normalized SU (M) generators 
Lk obeying ^ 1, ij, = L^, trL^ = 0, tr(ifc, L;) = M5ki 
|7|]. With \pk\ < 1 the hermitean matrix 

P=jj{'^ + PkLk) (13) 

obeys two of the requirements for a density matrix: tip — 1, 
trp2 = {P + l)/M < 1. The two-level observables A^''^ are 
represented as hermitean M x M matrices or "operators" 

i«=L, , (^W)=tr(piW)-p,. (14) 

Their expectation values can be computed according to the 
rule of quantum mechanics in terms of the density matrix. 
In particular, we consider the diagonal operators 

Li = dzag(l, 1, — 1, — 1) , L2 ~ diag{l, —1,1, —I), 
L3 = dzag(l, -1,-1,1). (15) 

The expectation values pi , P2 , Pa determine the diagonal 
elements of the density matrix 

1, X 1, N 

Pii = ^(1 + Pi+ P2+ Ps) , P22 = ^(1 + Pi - P2 - Ps), 

P33 = ^(1 - Pl+ P2- Ps) , Pii = ^(1 - Pi - P2 + Ps)- 

(16) 

While the purity constraint P < 3 assures paa < 1, the 
positivity condition paa > 0, which is necessary for a valid 
density matrix, is not obeyed automatically. (For exam- 
ple, pi = p2 = P3 = -1 ■ Pfc>4 = 0, obeys P = 3 but 
would imply pn = —1/2.) We will therefore extend the 
purity constraint and impose further conditions on the al- 
lowed values of pk (beyond the condition ptpk < 3) 
that will guarantee paa ^ 0. The explicit form of this 
positivity constraint will not be important for the present 
paper and we refer for a detailed discussion to [2*], [7]. In 
short, we require that the eigenvalues of the density ma- 
trix should all be positive or zero. For a pure state 
this implies Pk = fk, with fk parameterizing the homoge- 
neous space SU{M)/SU{M — 1) x C/(l), normalized with 
Y.khh = M-l. 
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4. Classical particle 

The specific properties of tfie "two-bit-particle" depend 
on the interpretation of the basis observables A^''\ For 
M — A, we could represent a classical particle by associ- 
ating the location L and the motion M to two commuting 
operators, for example 

L = Li= diag(l, 1, -1, -1) , M = L2= diag(l, -1, 1 - 1). 

(17) 

Thus L and M can have simultaneously sharp values. If 
no other observables are considered, only the values of pi 
and p2 are needed for the computation of (L) and (M) in 
the ensemble. A pure classical state has |pi| — \p2\ — 1- 
We may further consider an observable associated to the 
third diagonal generator L3 = L1L2 =diag(l, — 1, — 1, 1). 
It measures the correlation between motion and location, 
LM. This observable takes the value -1-1 for a particle 
located at r moving clockwise or a particle located at (l) 
moving anticlockwise, while LAI = — 1 holds if (r) and 
(Z) are exchanged. The most general statistical state of a 
two-bit classical particle is characterized by three numbers. 

pi = (L) , p2 = (M) , p3 = (LM). (18) 

For the example pi = P2 = , p3 = 1 , (LM) = 1, one 
does not know if the particle is located left or right or if 
it moves clockwise or anticlockwise. However, we know it 
moves "downwards" , either clockwise on the right or anti- 
clockwise on the left. As long as only the observables L, M 
and LM are used as particle observables and pi , P2 , Ps, de- 
scribe the state of the particle, we deal now with complete 
statistics, for which the measurement correlation {LM)m 
is expressed by the expectation value {LM) = p3, which is 
directly connected to the joint probabilities. The charac- 
teristic feature of a classical particle is that only diagonal 
operators are used for a description of the observables, and 
that the observables include both location and motion ob- 
servables. 

5. Quantum particle 

In contrast, for a description of a two-bit quantum par- 
ticle the characterization of the particle state also uses in- 
formation from observables corresponding to off-diagonal 
operators. We may now employ the three diagonal gen- 
erators Li,L2,L3 for a refined characterization of the lo- 
cation using four cells. (Imagine in our intuitive example 
that we place two further detectors at the location x = tt 
and X — —TT.) We may label the cells with a = 1 ... 4 
and "central locations" at x(a = 1) = 37r/4 , x{a — 2) = 
7r/4 , x{a = 3) = -7r/4 , x{a = 4) = -37r/4. The particle 
numbers iV*^"^ obey 



iv(i) = 




Li 


+ L2 


+ L3) 






Li 


-L2 


-L3) 


7V(3) = 




Li 


+ L2 


-L3) 






Li 


-L2 


+ L3) 



(19) 



and the location of the particle is given by X according to 
eq. dl]) , with eq. equivalent to eq. ^ . The observable 
Li takes positive values if the particle is on the right (x > 
0), and negative values for a particle on the left {x < 0), 
according to 

Li = 7VW-KiV(2)-iv(3)-7V(4). (20) 

Similarly, L3 is positive (negative) for a particle in the lower 
half, > n/2 (upper half, \x\ < n/2), 

Ls = N^^^ - N^^^ ~ N^^^ + (21) 

and L2 is positive or negative 

L2 = N'-^^ - iV(2) + iV(3) - Nii) (22) 

in one or the other of the two complementary diagonal 
regions of the circle. 

The off-diagonal generators L4 . . . L15 may be split into 
real generators (L^, = Lfe), namely L4 = 1 ® ti , Lq ~ 

T3 (g) Tl , Ls = n «) 1 , Lio = Tl (g) T3 , L12 = Tl (g) Tl , Ll4 = 

—T2®T2 and imaginary generators (i^ = —Lk), namely 

L5 = 1 <8) T2, L7 = T3 (g) T2 , Lg = T2 (g 1 , Ln = T2 g) 

T3 , Lis. = Ti g) T2 , L15 = T2 g) Ti. (In this language 
one has Li = r3 g) 1 , L2 = 1 ® T3 , L3 = r3 g) T3.) We 
may employ the six imaginary generators for generalized 
motion observables. One motion observable Mi can be 
associated with — £5, such that L5 = 1 for anticlockwise 
motion, and ^5 = — 1 for clockwise motion. We observe 
that L5 and Li commute such that there can be states 
where they have simultaneously sharp values. Indeed, a 
state with pi = (Li) = 1, ps = (L5) = —1 characterizes 
the particle as being in the right half of the circle and 
moving clockwise. In this special case we can recover a 
classical particle where a location and a motion observable 
are simultaneously measurable (even though with reduced 
precision, since only one bit can be used for the location 
measurement). Indeed, we may employ a change of basis 
for the generators Lk by unitary transformations (with a 
simultaneous transformation of p such that tr(pA) remains 
invariant). This can be used in order to transform the 
triplet of generators (Li, —L<^, —L^) into (Li, L2, L3), and 
we see the equivalence with the previous description of a 
classical particle if no information on further observables 
is used. This allows an interpretation of the observable 
associated to Lj. For Lt — 1 the particle is either on the 
right moving anticlockwise or on the left moving clockwise 
(the particle moves "upwards"), while for iy = — 1 it either 
moves clockwise on the right or anticlockwise on the left 
("downwards" motion). 

The observable Mi associated to Mi = — L5 is not the 
only motion observable. A different motion observable is 
associated to Al2 = -^13, with positive L13 describing clock- 
wise motion. This operator commutes with L3, such that 
we can now simultaneously answer the question if the par- 
ticle is in the lower half (L3 = 1) or in the upper half 
(L3 — —1), and if it is moving clockwise (L13 = 1) or 
anticlockwise L13 = —1. Again, the triplet of operators 
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{Ls, Li3, Li5 = L3L13) is equivalent by unitary transfor- 
mations to {Li, L2, L^) and we can describe a classical par- 
ticle. The third set of motion type observables, — Lg and 
— Lii = —L2Lq, commutes with L2- Its interpretation is 
less intuitive since it involves disconnected regions for the 
location in the ring. 

The two motion operators Mi and M2 commute (with 
M1M2 = —Lg). We may define an angular momentum 
observable 

M=-(M,+M.) = -(L,3-i5) = ^ V -1 I 1 

\ 1 0-10/ 
(23) 

This operator has three different eigenvalues, m = 0, ±1, 
as characteristic for angular momentum in quantum me- 
chanics in appropriate units. We observe that M does not 
commute with anyone of the location operators Li,L2,Li 
or with X. 

Pure states correspond to classical ensembles for which 
the density matrix obeys = p. This requires a maximal 
purity P = 3. Any pure state density matrix can be written 
as 

p = C/piC/t , pi = diag (1,0,0,0), (24) 

with U a suitable unitary matrix C/^J7 =1. For pure states 
we may introduce the quantum mechanical wave hmction 
as a normalized complex M-component vector ip, ijj'^tp = 1, 
according to 

= Uafilpl , i-lpn)a = <5na , P«/3 = Tpalpp, (25) 

such that the expectation value of an observable A can be 
computed using the associated operator A with the stan- 
dard quantum mechanical rule 

(A) = ip'fAip. (26) 

The eigcnstates of the location operator X are the eigen- 
states of the local particle numbers tpn, with eigenvalues 
Xn = {5 — 2n)7r/4. For all these "localized states" the 
average of the angular momentum observable vanishes, 
tp'<M^ = 0. 

The eigenstates of the angular momentum operator are 
given by 

1 r._ , , , (5 — 201)77 , , 

tPa = —exp[imx[a)\ , x{a) = , (27) 

with m = 0, ±1, ±2. We notice that rh = +2 and 
m = — 2 leads to the same state ip = ±|(1, — 1, 1, — 1) up 
to a minus sign. It corresponds to a zero eigenvalue of M. 
For |m| < 1 one has rh = m. For the angular momentum 
eigenstates one finds 

(AtM) = V)^iV("V =\, {X)= ^^X^p = 0. (28) 

It is amazing how many characteristic features of quan- 
tum mechanics are already visible for our simple "two bit 
particle" ! 



Indeed, wc can now characterize a two-bit quantum par- 
ticle by a subsystem with the following properties: (i) All 
particle observables have a spectrum of at most four dif- 
ferent values, corresponding to two yes/no questions, (ii) 
The state of the subsystem cannot be fully characterized 
by the expectation values for three "commuting observ- 
ables" , for which joint probabilities are available. It needs 
the specification of expectation values of further observ- 
ables. (iii) The statistical system describing the particle 
alone is incomplete. Joint probabilities are not available 
for all pairs of observables. (iv) The location and motion 
observables use a maximal resolution consistent with two 
bits. In consequence, they cannot be associated to com- 
muting operators. 

Finally, we may describe the classical particle within the 
quantum mechanical formalism as "classical particle eigen- 
states" which have simultaneously sharp values of Li and 
Ml (and therefore also of L^). The eigenstates of Li = +1 
obey = {tp 1, 11^2, 0,0). The simultaneous eigenstates with 
Ml = ±1 read ijj = (1, =Fi, 0, 0)/\/2. (For the eigenstates 
with Li = —1 one exchanges the upper two with the lower 
two components.) In a quantum mechanical language, the 
classical particle states are entangled states between the 
location eigenstates ijji and ■02 • 

6. Classical statistical ensemble for two-bit-particle 

There are many possible classical statistical ensembles 
that realize a two-bit-particle. They typically differ in the 
properties of the environment, while they lead to identical 
results for the subsystem which characterizes the two-bit- 
particle. As a first explicit example for a classical statistical 
ensemble we may consider 2^^ classical states r which are 
characctrizcd by ordered sequences {(7k} of fifteen discrete 
variables ak = ±1. We specify the state of the subsystem 
by the expectation values of the fifteen observables ak, 

For the probabilities Pr = p({ffe}) of the classical statisti- 
cal ensemble we choose 

P{Wk}) = Ps {Wk}) + Spe {Wk}) (30) 

with 

15 

PsiWk}) =2-'^l[{l + Pkak). (31) 
fe=i 

The part Spe obeys 

J2^lSPe{{^k})=0, ^(5pe(K})=0, (32) 
{(^k} {c^k} 

such that it only matters for the environment, without in- 
fluencing the state of the subsystem. This example can 
be easily generalized by extending the number of classi- 
cal states, using an additional index ( with t = ({cTfe}, C)- 
The probability Ps({crfe}) is multiplied by Ps{C) obeying 
JZrPsiO — 1> ^Pe depends on ( in addition to {crfe}- 
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The sums in eq. ([29)) . (p2)) became now sums over {crfe} and 
(. For example, the sequences of hits discussed in sect. IIIII 
can be described in this way. 

We may associate to each classical observable ak the op- 
erator Lfe and define the density matrix 

P=lil + PkLk). (33) 

This guarantees the quantum rule for the expectation val- 
ues 

(A) = tr(ip) (34) 

ii A — (Jk and A — Lk- We constrain the "purity" 
P ~ J2k PkPk < 3 and impose further on the allowed val- 
ues of pk the positivity constraint that all eigenvalues of p 
should be positive or vanish. Then p in eq. (1331) has all 
the properties of a density matrix for a four-state quantum 
system. A unitary time evolution of the density matrix p 
can be achieved by a suitable time evolution of the classical 
probability distribution (j30p . as realized by an appropriate 
rotation among the pk, while the evolution of dpe is arbi- 
trary as long as the constraints (|32|) are respected. 

One may realize the occupation numbers (jl9p in the 
classical statistical ensemble by appropriate sums of cti , <72 
and CT3. In general, these sums could take the values 
(1,^,0,—^), instead of the allowed values (1,0) for the 
occupation numbers. However, we may impose on the al- 
lowed classical probability distributions pr a further re- 
striction such that (Ti , (72 and (73 form a "comeasurable bit 
chain" This means that their expectation values obey 
the relation 

(cti • (T3) = (0-2) , (cti • (T2) = (0-3) , (0-2 ■ c^s) = (c^i), 

(cTi • CT2 • CT3) = 1- (35) 

Comeasurable bit chains are possible if the associated 
quantum operators commute. The relations p5p are then 
realized also on the operator level 

L1L3 — L2 , L1L2 = L3 , ^2-^3 — Li , L1L2L3 — 1. (36) 

The relations ([55)1 imply that the observables A^^\ to which 
the operators N'^^^ are associated, behave as projectors 

((A(^'))P) = (37) 

This may be verified explicitely for p = 2,3, and the use of 
the identity 

(A(^))4 = {A^^^f + i(yl(^'))2 - i(A(^')). (38) 

The relation ([57)1 implies that the probabilities for find- 
ing for the sum the values ±i vanishes, such that the 
spectrum of the observables A^^^ contains indeed only the 
values 1 and 0. Further comeasurable bit chains can be 
constructed for the motion observables (— (75, uis, — us), in 
complete analogy to (cti, CT2, 0-3). We can also implement 
a comeasurable bit chain (cri, — (75, — cry) for the classical 



particle. For a two-bit-particle or four-state quantum me- 
chanics the maximal numbers of members of a bit chain is 
three. 

Already on the level of two-bit subsystems we have con- 
structed classical statistical ensembles that can describe a 
quantum particle as well as a classical particle. The "na- 
ture" of the particle depends on what observables are used 
for the description. In our 2-bit example, the angular mo- 
mentum operator for a quantum particle is M (eq. (j23p ). 
while for a classical particle one uses Mi = — L5. For a 
classical particle, typically not all of the information con- 
tained in the density matrix p is used - with the exception 
of states where = 0, fc 7^ 1, 5, 7. While the location in- 
fomation of a quantum particle uses the three expectation 
values pi,p2,P3j only pi is used for a classical particle. In 
the next sections we will generalize these findings by con- 
sidering a continuum limit where the number of bits goes 
to infinity, thereby increasing the resolution of the location 
observable X and the range of the angular momentum M . 

IV. PARTICLE MOTION 

In this section we extend the discussion from the "one- 
bit-particle" and "two-bit-particle", where all observables 
have a discrete spectrum, to quantum and classical parti- 
cles described by continuous position and momentum ob- 
servables. For the quantum particle we will obtain Heisen- 
berg's uncertainty relation and non-commuting position 
and momentum operators. In contrast, for a classical parti- 
cle position and momentum can both be measured sharply 
and the associated operators commute. 

1. Continuous observables 

Let us next increase the number of bits Q and the cor- 
responding number of quantum states M — 2^ . We finally 
will be interested in the continuum limit M ^ 00. We now 
consider a basis of — 1 two-level observables A'^^^ for 
describing the state of the particle by eq. (|lip . Equations 
(1121114^ have been already formulated for general M (now 
P < M — 1) and we keep the same normalization for the 
SU{M) generators L^. Our quantum mechanical formal- 
ism holds for general M. For the particle numbers TV*^"^ 
we use eq. (O, with a = 1 ... A/. Defining 

x{a) = {M + l-2a)^^ (39) 

we infer the "quantum location operator" X as 

X = ^x{a)N^"\ (40) 

a 

Thus X remains an angular variable on the circle, — tt < a; < 
TT, but the resolution increases due to the higher number 
of possible eigenstates. 

For M — > 00 the location becomes a continuous variable. 
We may switch to a continuum notation where V'q is re- 
placed by the continuous wave function ^l){x). The density 
matrix becomes a function of two coordinates, p{x, y). For 
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the special case of a pure state it reads p{x, y) = ip{x)ip* (y)- 
Similarly, the operators become matrices in position space, 
A{x,y). In particular, the particle number at location z 
reads 

7V(")(x,y) =(5(x-z)(5(y-z). (41) 
The quantum position operator takes the form 

X{x,y)^ j dz zN{z){x,y)=x5{x-y), (42) 

with expectation value 



(X) = / / X{x,y)p{y,x) = / xp{x,x). (43) 

J X J y <y X 

We can interprete p{x,x) as the probability w{x) to find 
the particle at location x. It obeys < wix) < 1 and 
the condition trp = 1 assures the proper normalization 
w{x) = 1. For a pure state we have the usual rule 



{X) = / '4'* {x)xip{x) , w{x) — Tp* {x)ip{x). 



(44) 



2. Momentum and angular momentum for 
quantum particle 

Consider next the angular momentum operator M for a 
general number of quantum states M. (There should be 
no confusion between the similar symbols.) It reads 



M 



M, 



(45) 



7 = 1 



Here 7 is considered as an index modulo Af, i.e. 7 = + 1 
corresponds to 7 = 1. The normalization Af depends on 
the units for M. We choose 



.r ^ . 27r 

Af ^2 sin — 

M 



such that 



exp ix{a) 



(46) 



(47) 



is an eigenstate of M with eigenvalue m — 1. One finds for 
the eigenvalue rh{m) of M, with eigenfunctions 

V'i™^ = ^ exp [imxia)] , - y < m < y (48) 



the spectrum 



^ , , sin(27rm/M) 

rh(m) = ; — — . 

^ ' sin(27r/M) 



(49) 



For M ^ 00 the small values of rh become the usual 
integer values of angular momentum, i7i{m) — > m. For 
m — Af/4 the function rh{Tri) reaches its maximum m,iiax = 
1/ sin(27r/A'/) « Af/27r, and then decreases towards zero 
for TO — > Af/2. We see that the range of possible angular 



momenta increases ~ M, in correspondence with the in- 
creased angular resolution. Ax ~ 1/M. The appearance of 
small eigenvalues to for \m\ Af/2 is an artefact of the lat- 
tice formulation (similar to the fermion doubling in lattice 
gauge theories). They may be removed towards large val- 
ues by the use of an improved angular momentum operator. 
We will assume that this is done and consider for the con- 
tinuum limit only smooth functions with \m\ <^ M/A. In 
the continuum limit we find integer eigenvalues of angular 
momentum. The continuum version of M is a differential 
operator 



M{x,y) = -iS{x-y) 



d_ 
dy' 



(50) 



If we further restore standard units by multiplying M with 
h (so far we have used h = 1) one recovers the well known 
commutation relation between an angular position opera- 
tor X and the angular momentum operator M 



[X,M] ^ihS{x-y). 



(51) 



The continuum limit can also be taken in a different way. 
While increasing the resolution we may simultaneously in- 
crease the "volume" of the circle (i.e. the circumference) 
to infinity, and consider only a finite region in this infinite 
volume. This essentially amounts to a change of units and 
to a restriction of observables to a shrinking region in x. 
We may choose some length unit and perform a rescaling 
x{a) — 2t:x' /I, such that the new location variable extends 
from —1/2 < x' < 1/2, with I = Ma the volume and a the 
lattice distance. With X' = {1/2'k)X , P = {2TTh/l)M, the 
new variables X, M' read 

X' = x'Six' -y') , P = -mix' - y') — . (52) 

dy' 

We interprete now P as momentum - in one dimension 
there is no difference between momentum in a periodic vol- 
ume and angular momentum (up to units). The relation 
([5T|) becomes the well known Heisenberg's uncertainty re- 
lation between position and momentum operators (we drop 
the primes oti X ,x from now on) 



[X,P] = m{x-y). 



(53) 



For a finite I the momenta remain discrete, with p = 
2'k/I ( from here on we use again ?i = 1). A scaling 
I ~ ^/M,a ^ 1/-\/M, however, leads for Af — cx) to con- 
tinuous momentum and location variables. Restricting the 
discussion to a fixed finite range Ax, with Ax /I — for 
Z — > 00, the periodicity of x can be neglected. We end with 
the standard setting for a particle in quantum mechanics. 
In this version the generalization to three dimensions is 
straightforward. 

3. Commuting position and momentum for 
classical particle 

What happens for Af — > 00 with the "classical parti- 
cle state" that we have found in our quantum mechanical 
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formalism for M = 4? It is indeed straightforward to gen- 
eralize for large M the notion of a classical particle where 
both location and momentum can be measured. For this 
purpose we use half of the yes/no questions for the location, 
and the other half for the motion. (We assume here Q/2 
to be integer.) We can construct a coarse grained location 
operator as a %/M x ^/M matrix (VM = 2'^/'^). It 

follows the same rules as the previous quantum location op- 
erator Xm, but with less resolution since only Q/2 bits can 
be used. Similarly, we construct a coarse grained momen- 
tum operator P The classical location and momentum 

operators Xci and Pd are given by M x M matrices 

XrJ = ® 1^ , Pa - l^M ® P^. (54) 

Obviously, classical location and momentum commute 



^d,-Pd] =0. 



(55) 



This construction generalizes the classical location and 
motion for M — A. Indeed X2 — T3 and P2 = — T2 repro- 
duce Xci = T3 I = Li , Fci = 1 ® (-T2) = -^5- In 
the limit M — > cxd both Xd and Pd become continuous 
operators. This demonstrates clearly that our setting of 
describing a particle as a subsystem of a classical statisti- 
cal ensemble can account for both quantum and classical 
particles! 

There are many different ways of taking the limit M — > 
00 which result in the classical commutation relation (|55p , 
and similar for the quantum commutator (1531) . At the end, 
only the commutation relation matters for the distinction 
between quantum and classical particles. As an example, 
we may represent the location observable for both a clas- 
sical and a quantum particle in the \/M-dimensional sub- 
space 



Xq = Xd = X^ ® 1^ 



'M- 



(56) 



The momentum operators can then be represented as act- 
ing in different subspaces 



Pq = K 



M 



1 , Pd = 1 «) p, 



M- 



(57) 



In other words, the M-dimensional position space can be 
parameterized by a pair (11,0:2), with representations 



Xr 



X,. 



^1 , Pq 



-ih- 



d_ 

dxi 



-ih 



d_ 

"0X2 



(58) 



In this representation the coordinate X2 becomes irrelevant 
for the observable describing a quantum particle. 



PARTICLE CORRELATION FOR QUANTUM 
PARTICLE 



So far we have only discussed the diagonal and the imag- 
inary part of the density matrix for quantum particles. In 
our basis they are related in a general sense to the location 
and motion of the particle. The real off-diagonal part of 
the operators and the associated part of the density matrix 



describes a type of single particle correlation. It general- 
izes the observable C for the one bit system which takes 
the value -1-1 if the particle has a tendency to be both right 
and left (correlation), and —1 if the tendency is towards 
mutual exclusion of the two locations. Similarly, for the 
two bit system a density matrix 



1, 1, 0, 

1 [ 1, 1, 0, 

2 I 0, 0, 0, 

0, 0, 0, 0. 



(59) 



indicates a tendency that the particle is simultaneously in 
the first quarter of the circle (x — 37r/4) and in the sec- 
ond {x = 7r/4). One finds {M^^) = (M^^) = 1/2, while 
(iV^'^)) = (N'^^'I) = 0. However, the state is not a mixed 
state with equal probability of the particle to be in one 
of the first two quarters. Such a mixed state would cor- 
respond to a density matrix p = diag{l/2, 1/2,0,0), with 
trp^ — 1/2, whereas the correlated state ([5^ is a pure 
state with — p. Indeed, the density matrix ([M)) can be 
associated to a wave function tp — {tpi + 'tp2)/V2- 

1. Extended wave functions as eigenstates of 
correlation operators 

This can be generalized immediately to the continuum 
limit with infinitely many bits. Consider a wave function 
which corresponds to a Gaussian wave packet 

^(x) = i2nA^)-i exp j-^^^j • (60) 



exp j 

The corresponding density matrix 
p{x,y) — (27rA^)~5exp 
{x - yf 



(x + y- 2xf 
8A2 



X exp 



8A2 



(61) 



shows a nonvanishing correlation between distant x and y. 
For example, the operator A = 5{x — y — 2a) can be used 
in order to test the shape of the disctribution. The density 
matrix (j6ip leads to a nonvanishing expectation value of 
the associated observable 

„2 



{A) = I P{x, y)S{x - 2/ - 2a) = exp - 



xy 



2A2 



(62) 



which vanishes for large a with a characteristic width given 
by A 

The expectation values of the location observable X and 
its dispersion read 

{X)^x , {{X~x))^A\ (63) 

Thus eq. (|61l) describes a situation where physics in neigh- 
boring regions, \x — y\/2 < A, is strongly correlated. We 
may associate Re[p{x,y)) with a correlation function. 

The Gaussian wave function ((60|) is an eigenstate of the 
family of hermitean correlation operators 

CA{x;x,y)^-\^ (64) 



47rA2 



1 



X exp ■ 



4A2 



a{x - yf + 



y 
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obeying 



with 



dy Ca{x; x, y)'ipiy) = CAipix), 



(65) 



4A4i2 



(72) 



with eigenvalues 



CA 



1 + 2a Vtt/ 



a\i/2 1 
A' 



(66) 



We have normahzed Ca such that in the Umit A — it 
becomes the occupation number at location x 

hin Ca = = S{x ~ x)6{y - x). (67) 

However, for any A > the operators C'a are real symmet- 
ric operators of the particle correlation type. The family 
of operators depends on a parameter a. For a — >■ cx) the 
correlation operator Ca becomes proportional to the unit 
operator and independent of x, while for a — ^ it becomes 
proportional to J-distribution for the center of mass coor- 
dinate z = (x + y)/2, i.e. ~ 6{z — x), while the relative 
coordinate plays no role. 

2. Wave packets 

As well known from quantum mechanics, a moving quan- 
tum particle can be described by a wave packet. We reca- 
pitulate here the basic properties in order to facilitate the 
comparison with wave packets for classical particles that 
will be discussed in the next section. Let us consider a 
wave function 



(68) 



where we choose for definiteness a nonrelativistic free par- 
ticle with w(p) = /2m, and a Gaussian wave packet, 



In the limit Ap — > one recovers the plane wave density 
matrix (eigenstate of momentum) 



(73) 



Pp{x, y) \ -Ap exp[ip(x - y)]. 



while for Ap — > oo one finds the limit of a sharply located 
particle 

-^^(^ (74) 



Using the momentum operator (j52p . P{y,x) — ^i5{y 



x) -§- , the density matrix (1711) describes a time independent 



momentum distribution 

{F{p)) = / {F{P)){y,x)p{x,y-t) 



(75) 



= lim / F{~idx)p{x,y;t) 



dp 



F{p)p{p,p\ t) 



dp 



F{p)A\p). 



27r ' J 27r 

Here we employ the Fourier transform of the density matrix 



p[p, q; t) 



pix,y;t)e-^P^e^'^y 



A{p)Aiq)e^p{^iq^-p')}, (76) 



and {P)^p, ((F-p)2)=A2. 

On the other hand, functions of the location observable 
obey 



In Fourier space this yields 



p 



(69) 



{F{X)) = I F{x)p{x,x-t) 
Pi 



(77) 



Ap / F{x) exp <^ -2A^ 



pt 



ip{p,t) = A(p)exp -i 



.p^t 
'2m 



(70) 



Again, the wave packet ()68p describes a state with non- 
vanishing particle correlation for y ^ x, similar to the static 
Gaussian distribution ([60t . 

The corresponding pure state density matrix reads 



p{x,y;t) 



-Apexp<^ -2A^ 



X + y pt 



X exp 



A2 



(71) 



X exp <.i{x — y) 



P + 



AAmt(x + y pt 



This implies for the location and its dispersion 

{X)^x{t)^^, {{X^x{t)f) = ^. 

rn 4A„ 



(78) 



For i — > cx) the vanishing of Ap{t oo) implies a 
diverging dispersion, with Ap — > m/{2Apt) and 



p{x,y;t-^ oo) 
m 



InApt 



exp ■ 



2A2 



m{x + y) 
2t 



X exp ■ 



1 / m{x — y) 



2AI 



2t 



mlx + y). , 
X exp<;i (x~y) 



(79) 
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We finally display the Wigner representation of the 
density matrix ([7T|) (Wigner quasi-probability distribution) 
by using the center of mass coordinate z ^ {x + y)/2 and 
performing a Fourier transformation with respect to the 
relative coordinate x — w, 



Pw{z,q;t) 



= 2exp|-2Af('^-^)'>| 
xexp{-2A^(.-£yj. 



(80) 



For this particular case of a free particle Pw{z,q) consti- 
tutes a probability density in the classical phase space of 
location and momentum. It is real, positive, and normal- 
ized, 



dz 



dq 
2^ 



Pw{z,q;t) = 1. 



(81) 



Marginalizing over the center of mass distribution yields 
the normalized probability distribution of momenta 



p{q) = / dzp^{z,q) 



/27r 



■ exp 



1 



2A2' 



iq~pr}=A'iq), 
... (^2) 

while integrating over momenta we find the distribution of 
the center of mass coordinate around the average classical 
trajectory z = pt, 

dq 



P{z) = / -^Pw{z,q) 



Ap exp ■ 



-2a; I z 



m 



(83) 



We conclude that a quantum particle is typically de- 
scribed by nontrivial particle correlations. In the Wigner 
representation the density matrix for our particular wave 
packet for the free particle defines a probability distribu- 
tion in classical phase space of location and momentum. It 
is well known that this property does not hold for a general 
density matrix of a quantum particle - the Wigner function 
Pw can be negative in certain regions of phase space. Fur- 
thermore, the positivity of p^ may not be preserved by the 
time evolution for an interacting particle. 



VI. QUANTUM FORMALISM FOR CLASSICAL 
PARTICLE 

Let us now turn again to the possibility of realizing a 
probability distribution for a "particle" with arbitrarily ac- 
curate location and momentum. Since we can associate 
a quantum mechanical density matrix to this probability 
distribution, this constitutes a way to describe a "classical 
particle" , with sharp location and momentum, within the 
formalism of quantum mechanics. Of course, the relevant 



location and momentum operators must now commute, cf. 
eq. (1551) . We emphasize that this classical particle is not 
the usual classical limit of quantum mechanics. In this 
section we rather develop the formalism of quantum me- 
chanics for a microscopic classical particle. 

This part does actually not need our previous discussion 
how quantum mechanics arises from classical statistics. It 
constitutes a self-consistent probabilistic description of a 
classical particle. The basic concepts are the probability 
distribution in phase space and its dynamics. The latter is 
formulated as a time evolution equation for the probability 
distribution. This replaces the notion of trajectories and 
Newton's laws as basic concepts for classical particles. We 
first discuss the Hamiltonian which leads to a time evo- 
lution which corresponds to classical trajectories for par- 
ticles. In the presence of a potential, the time evolution 
according to the Schrodinger equation associated to this 
Hamiltonian changes the shape of a wave packet. We es- 
tablish that the probability density in phase space, which 
obtains from a given initial wave function, describes pre- 
cisely the classical probability distribution evolving accord- 
ing to the Liouville equation. In the traditional approach 
the latter follows if non-interacting particles move on clas- 
sical trajectories in a potential, with a distribution of initial 
conditions given by the squared initial wave function. 

A new basic concept in our treatment will be the wave 
function for a classical particle, which shares many formal 
features with the wave function for a quantum particle. Im- 
portant distinctions remain, however. The classical wave 
function depends on both position and momentum and is 
a real function. 

1. Wave function and quantum Hamiltonian for 
classical particles 

The direct product structure ([5^ for classical location 
and momentum operators implies that for each location 
X we can assign further quantum numbers. Due to the 
vanishing of the commutator (j55p we can choose the wave 
function to be simultaneously an eigenfunction of Xd and 
Pel- A general wave function depends then on both the 
variables x and p, with 

Xciip{x,p,t) ^ xip{x,p,t) , Pci'ip{x,p,t) ^ pip{x,p,t). 

(84) 

The structure of quantum mechanics perfectly allows for 
commuting position and momentum operators, which are 
represented here as 

Xci = xS{x - y) , Pel = pS{p - q), 

[X,i,P,i]^0. (85) 

(In our normalization S{p — q) stands in d dimensions for 
{2tt)'^S'^{p — q).) For the commuting operators Xci and Pd 
we can compute expectation values of composite observ- 
ables 

{F(Xd.Pci))^ j mx,p)\^F{x,p). (86) 

We can therefore identify the classical probability density 
in phase space, w(x,p), as 

w{x,p) = \^{x,p)\\ (87) 
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We will see that ■0 can be taken real, w = tp'^, such that 
■0 = i-x/w is given by w up to a sign. 

The wave function ip{x,p) is the central object for casting 
a probabilistic theory of classical particles into the quan- 
tum formalism. It obeys the generalized Schrodinger equa- 
tion 



idtip{x,p) = H->Jj{x,p), 



with H the "quantum Hamiltonian" . For an appropriate 
choice of H the time evolution of the probability distri- 
bution w — \ip\'^ is determined by the Liouville equation, 
as standard for classical particles. All the usual quantum 
rules for the computation of expectation values of observ- 
ables apply. 

The difference between a classical and a quantum par- 
ticle does not arise from the different formal structure be- 
tween quantum mechanics and classical statistics. We have 
seen before how to implement a quantum particle within 
the conceptual framework of classical statistics, and we es- 
tablish now how to describe a classical particle within the 
formalism of quantum mechanics. The difference between 
quantum or classical behavior of a particle is rather due 
to the different dynamics or, in other words, to different 
Hamiltonians which use different types of operators. 

As a special limiting case we want to describe within 
the formalism of quantum mechanics a particle with sharp 
location and momentum which follows a classical trajec- 
tory. For this purpose we cannot use a Hamiltonian 
Hci — V{Xci) + P^i/2m, since in this case both Xci and 
Pel would commute with Hd , resulting in conserved x and 
p. We may, however, use a Hamiltonian which involves 
both Pel and a derivative operator which corresponds to 
the usual momentum operator in quantum mechanics. In 
this section we concentrate on a Hamiltonian which reads 
in the (x,p) representation 



ih d ,dV d 

Hl = Ptt- + *"^Tr, 

m ox ox op 



(89) 



with V — V{Xci) represented as V{x) and omitting the 
(5-functions. Other choices of the Hamiltonian are possible, 
as the one leading to the time evolution for a quantum par- 
ticle which will be discussed in sect. IVIIH or the extended 
Hamiltonian investigated in appendix A. 

The choice of the Hamiltonian Hl leads to the commu- 
tation relations 



[HL,Xe 



ih 



Pel , [HL,Pci]^ih 



dV 
dUi' 



(90) 



(We recall that Hl denotes here the Hamiltonian operator 
in the quantum formalism and should not be confounded 
with the classical Hamiltonian in classical mechanics.) A 
standard quantum mechanical calculation yields then the 
time evolution for the expectation values Xd = {Xd} and 
pel = {Pel), according to 



OtXd = , OtPel = -{-rr-)- 

m ox 



(91) 



These are the same equations for the expectation values as 
for a quantum particle. They reduce to classical trajecto- 
ries if {dV/dx) can be replaced by {dV/dx){xd)- This can 
be realized by a proability distribution which is sharp in 
position space. We emphasize that the second eq. (|9T|) is 
a perfectly classical statistical equation if we describe the 
time evolution for a distribution of initial conditions. 

2. Time evolution of free classical wave packet 

We will see below that the Hamiltonian Hl leads to the 
Liouville equation for the time evolution of the probability 
density w{x,p). It is useful, however, to understand the 
time evolution of the wave function ip^x^p), as arising from 
the standard quantum formalism. Let us first consider the 
case of a free particle, dV/dx — 0. We start with an initial 
wave function 



(x — x^'^ 

ip{x,p) = {A^Ap)-i cxp\ ("^^P 



which is normalized according to 



{p - pf 



4A2 



(92) 
(93) 



In the limit A^^ — > 0, Ap this is a simultaneous eigen- 
state of location and momentum, with eigenvalues x and 
p. The time derivative according to the Hamiltonian ([5^ 
reads 



dtip = -iHlV = T^, 

m ox 



(94) 



where we use now again h — 1. Inserting a wave function 
of the type dH]) yields 



dtip 



p{x — x) 
2mA2 



(95) 



On the other hand, if we assume constant A2;,Aj,, with 
only X depending on time, we find 



1 



dti^ = ^^(^ ~ x){dtx)ip. 



(96) 



The expressions ([M]) and ([M)) coincide if we choose x de- 
pending on p and t, x{p,t) — xq + pt/m. For — > 
, Ap — > we can indeed describe a free particle 
with sharp location and momentum, moving on a clas- 
sical trajectory. As mentioned before, we can interprete 
Tp* {x,p)ip{x,p) as a probability for the particle in phase 
space {x,p). Comparing the wave function (j92p with the 
density matrix for a quantum particle in the Wigner repre- 
sentation ([50]) . we find that ip*tp coincides with eq. ((5D|) for 
Ax = l/(2Ap). In contrast to the quantum particle, how- 
ever, we can now choose both Aa; and Ap to take arbitrary 
values - they are not restricted by Heisenberg's uncertainty 
relation. 

One may try to generalize the free particle construc- 
tion in the presence of a potential, by using in eq. (I92p 
for p a location and time dependent mean value obeying 
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dtp{x,t) = —dV/dx. This will, however, not yield a solu- 
tion of the Schrodinger equation 



P 



dxi^ + idxV)dp'4>, 



due to the additional contributions involving 



(97) 



(98) 



For constant Ax,Ap in eq. (1921) the l.h.s. of the 
Schrodinger equation yields 

fl; P{x-x) p-pdV 

dt^^ ^T^dtx + ^^dtp^ ^ -^^-g^^ (99) 



2A2 



while one finds for the r.h.s. 



2mAi 



p dV 

— dxip + -^dpil) = 

m ox 



p{x — x) p — p dV 



2mA2 



2A2 



dx 



(100) 



with 



R 



pjp - P) 
2mAj 



dxP- 



2AI 



-dxVdpX 



p{p — p) d'^V 
dx"^ 



dV 



2toA2 



2mA2 dx 



t. 



(101) 



We conclude that a wave packet given by eq. where 
the time evolution arises only from x{x,p,t) and p{x,p,t), 
obeys the Schrodinger equation only for R = 0. In the 
presence of a potential this will, in general, not be the case. 
The time evolution of the classical wave function according 
to the Schrodinger equation changes the shape of ■0 beyond 
the form (|92]). 



3. Classical phase space distribution from 
Schrodinger equation 

For a classical probability distribution for pointlike parti- 
cles following classical trajectories in a potential, an initial 
Gaussian distribution will not remain Gaussian. We may 
start at time = with a real '0(x,p, 0) given by an ini- 
tial Gaussian distribution. Solving the Schrodinger equa- 
tion, the classical probability distribution w{x,p) at some 
later time obtains from the "wave function" ^ as given by 



il)(x,p,t) = {AxAp) ^exp' 



{xoix,p,t) ~ xY 



X exp • 



{po{x,p,t) -p) 
4A2 



(102) 



Here we invert the solution for the classical trajecto- 
ries for given initial values Xq^pq at t — 0, namely 
x{xo,po,t),p{xQ,pQ,t), in order to "extrapolate back" the 
initial values xo,po which correspond to given values of 
x{t),p{t). This defines xo{x,p,t) and po{x,p,t). We as- 
sume for simplicity that two different initial conditions re- 
sult in two different points in phase space {x{t),p{t)) for all 



t, such that the phase space trajectory is invertible every- 
where. For a real initial ^^:{xo,po) we observe that ipix,p,t) 
remains real and does not change sign during the time evo- 
lution (|102p . This may be different if we choose a Hamil- 
tonian different from H^. 

We have to show that the wave function (|102p is a so- 
lution of the Schrodinger equation with Hamiltonian(| 
For the wave function (I102p one obtains 



dt^j = - 



xo{x,p,t) - X 
Pt){x,p,t) -p 



dtxo, 



2AI 



dtPo,, 



V'- 



On the other hand, evaluating the Hamiltonian 
wave function (|102l) yields 



(103) 
for the 



iHLil^ = dxijj + dxVdpip 



m 



p{xo{x,p,t) - x) 



dxXOu 



p{po{x,P,t) -p) 



2mAl 

dxV{xo{x,p,t) - x) 
2Al 

dxV{pQ{x,p,t) -p) 

2A2 

p 



2mA2 




(104) 



We next use the observation that in the absence of an ex- 
plicit time dependence of the classical Hamiltonian a given 
initial value (a;o,po) can be connected to the trajectory at 
two different times ti and ^2 

{xQ,Po){ti,Xi,Pi) = {xo,Po){t2,X2,P2), (105) 

where 

Xl = X{ti) , X2 = X{t2) , Pi = p{ti) , P2 = p{t2) (106) 

obey for infinitesimal t2—ti the classical evolution equation 
for the trajectory (at fixed (a;o,po)) 



X2 — Xl 



P2 -Pi 



dx T) 

-j^ih-h) = -{t2-ti), 

at m 

dv dV 
^(t.-0--^(^.-M. (107) 



Expanding the r.h.s of eq. ()105|) around xi,pi one finds 
^tXo^^^^ 



dx dp 

dt ^ I"'' dt 



dtPou 



-—dxXf): J + dxVdpXo.^, 
dx dp 

-—dxPo\pt + dxVdpPo.f 



(108) 



Inserting these relations into eq. (jl03p we find eq. 
(|104l) . Thus the wave function (I102p . with XQ{x,p,t) and 
Po{x,p,t) determined by the classical trajectories, obeys 
the Schrodinger equation for the Hamiltonian 
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This finding generalizes for an arbitrary wave function 
of the form 

i){x,p,t) = i:(xQ{x,p,t),pQ{x,p,t)). (109) 
Indeed, the relation 



dxn 



dtXQ. 



dip 
dpo 



dtPo\x,p 



(110) 



shows that ip obeys the Schrodinger equation if 
XQ{x,p,t) , po{x,p,t) are determined by the classical tra- 
jectories according to eq. (jlOSp . using 

^^'^i'^'' " ' ^ ^^^^^ 

In consequence, the quantum mechanical probability 
{tjj*tjj){t,x,p) for finding a particle at time t at the loca- 
tion x with momentum p is exactly the same as the classi- 
cal probability for pointlike particles that follow a classical 
trajectory each. This holds provided that at some time 
to the initial classical probability distribution for the po- 
sitions of the classical particles in phase space is given by 
{tp*ip){tQ,xo,po). 

In summary, we have established a quantum mechan- 
ical description for classical particles. Using a Hilbert 
space with commuting position and momentum operators 
Xci , Pel , and the Hamiltonian (|89p , we find that states with 
simultaneously sharp location and momentum are possible. 
Furthermore, the probability distribution in phase space 
evolves exactly as for classical particles following classical 
trajectories. As long as we concentrate on "commuting ob- 
servables" , which can be written as functions of x and p, as 
the energy E = V{x) +p^ /2m, there is a one to one corre- 
spondence between our quantum mechanical model and the 
usual classical description of the evolution of distributions 
in phase space. Indeed, since L — iHl is a real operator 
containing only first derivatives, we can directly translate 
the Schrodinger equation for ip to s. similar time evolution 
equation for the probability density in phase space 



dtw{x,p) — —iHLw{x,p) = —Lw{x,p). 
This is the Liouville equation. 



(112) 



4. Quantum observables for classical particles 

The general solution of the Schrodinger equation for ip 
is given by 

ip{x,p,t) = tl;o{xQ{x,p,t),po{x,p,t)), (113) 

with 4'o{xoiPo) the initial wave function at some time tg, 
and (^xo{x,p,t),po{x,p,t) determined by "following back" 
the classical trajectories, as discussed above. In general, 
-0 could be a complex normalized vector in Hilbert space. 
Since the evolution equation dt^p — —Lip involves only real 
operators, the real and imaginary part of ip do not get 
mixed during the evolution. For a real initial wave function 



iPq{x,p) the wave function ip{x,p,t) remains real for all t. 
In general, we may write ip in terms of the probability 
density w{x,p,t) and a phase a{x,p,t), 

ip{x,p, t) — w^^^{x,p, t) exp{ia(x,p, t)}. (114) 

Only w^^^{x,p,t) is needed for the computation of ex- 
pectation values and correlations of arbitrary observables 
of the type A{x,p). Thus for all questions in classical 
physics the phase a is redundant. It only appears in the 
off-diagonal elements of the density matrix p{x,y,p,q) = 
ijj{x,p)ip* {y,q), while it drops out from the diagonal ele- 
ments p{x,x,p,p) — w{x,p). Only the diagonal elements 
of p influence the expectation values of diagonal operators 
A = A{x,p)d{x — y)S{p — q). For pure states we find for 
the modulus of the elements of the density matrix 



\p{x,P,y,q)\'^ = w{x,p)w{y,q). 



(115) 



Again the phase a plays no role. 

This changes if we ask further "quantum mechanical 
questions" that do not appear in classical physics. For 
example, is a hermitean operator and we may want 
to compute its expectation value. For the general ansatz 
(|114p it depends on the phase a{x,p) 



(Hr) = 



x,p 



— OxOL - -T^Opa 
m ox 



(116) 



and differs from the average energy (E) = {p^ /2ra + V{x)) . 
The additional information contained in the phase of the 
wave function would then be needed for the computation of 
expectation values of off-diagonal operators. We will, how- 
ever, mainly restrict the discussion to a real wave function 
ip{x,p). Then the phase can only take two values and tt, 
according to the sign of ijj. Instead of using w and a it is 
more convenient to use directly the real function ip. For 
real ip one finds (Hl) = 0. 

We can generalize our discussion of a "real time evolu- 
tion" for an arbitrary antisymmetric and purely imaginary 
Hamiltonian. One can always define "classical subsystems" 
by restricting the observables to a set which corresponds 
to mutually commuting operators. In a basis where these 
operators are diagonal the phase of the wave function be- 
comes irrelevant. (Note that the phase depends on quan- 
tum numbers as {x,p), such that we deal with different 
phases for each {x,p). This should not be confounded with 
an overall global phase of the wave function, which is al- 
ways irrelevant in quantum mechanics.) Under this angle 
the typical quantum mechanical features are connected to 
the off-diagonal operators. Their expectation values may 
depend on the additional phase information. Still, a real 
time evolution always allows a setting with real ip, such 
that the phase information concerns only the sign of ip. 
The situation is different for Hamiltonians which contain a 
part which is a real function or differential operator. We 
discuss an example of such an "extended Hamiltonian" in 
app. A. In the following we restrict the discussion to a real 
wave function ip{x,p) and purely imaginary H. 

Let us next consider the hermitean operators 



R = -idx. , ^.s = idr, 



(117) 
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which obey 



Keeping e, I and as parameters one finds 



[Ps,Xci] = -i , [Xs , Pel] = i , [Xs,Ps] - 0, 
[Ps,Pci] = , [X,,X,i]=0. 



(118) 



We can write the Hamiltonian in terms of these operators 



1 



Hl = -PciPs + V'{Xci)Xs 



(119) 



Both Xg and Pg are off-diagonal operators (in a basis where 
Xci and Pel are diagonal) and belong to the quantum ob- 
servables. For a general complex wave function given by 
eq. (jll4l) their expectation values depend on the phase a. 
We will see how for a real wave function ip{x,p) the expec- 
tation values can actually be expressed only in terms of the 
probability density w{x,p). 

The expectation values of Xs and Pg vanish by the ab- 
sence of boundary terms 



(Pa) = ^ix,p)dxip{x,p) = I d^w{x,p)^0, 

Jx,p ^ 



{Xg 



0. 



(120) 



Nevertheless, the squared operators have positive, in gen- 
eral nonzero, expectation values 

(P2) = _ /■ ^{x,p)dl^P{x,p) = / {d,ij{x,p)Y 
Jx.p Jx,p 



X.p 'J x,p 



(X-^) = 1 I {dp^x,p)y = - w{d,lnw)\il21) 

Here we use dxW = 2tpdx'i{^ , {dxw)"^ — 4:w{dxip)'^ and 
consider classical probability distributions w for which the 
r.h.s. of eq. (|12ip is well defined. For products 



vn pn 



(122) 



one finds {G„m) = for odd n or m, while for both n and m 
even we can express (Gnm) in terms of dx Inw and dp Inw 
without an explicit dependence on the sign of ^, similar 
to eq. (|12ip . We conclude that Xg and Pg are quantum 
observables for which operators of the type G can be com- 
puted in terms of the classical probability distribution w. 

5. Statistical observables 

The expectation value (Pg) in eq. (|121l) is a measure for 
the roughness of the probability distribution in position 
space. Consider the example 



j{x) = J w{x,p) = cos^ (^—^ exp 



2A2 



(123) 



with e chosen such that 



W{e^,l,A^) = / w{x) = 1. (124) 

J X 



1/2 



and we can easily compute 



ID 



= w- 



2A^ 



d 



5A2 

{XD = 2/X^{d\nW/dA?) 



1 + exp 



2A^ 



(125) 
(126) 



1 



4A2 exp(-2AV;2) 
l2~i + exp(-2A2/P 



For P //S? — >• one finds that {X'^i) A^ becomes inde- 
pendent of Z up to exponentially small corrections, while for 
Z^/A^ oo the dispersion reaches A^ up to power correc- 
tions, {X^i) — >• A^ — A^'/Z^. The situation is qualitatively 
similar for higher powers {X'^l''). 
On the other hand, we have 

{Ps) - j^{dxH^))\ 
i{j{x) = / ?/;(x,p) = ecos (^yj 



exp 



4A2 



(127) 



where we assume here for simplicity a factorized ^{x,p) 
t{p)ijj{x) with £2 = J^e'^{p). We find that 



(Pf 



Z2 l + exp(-2A2/Z2) 4A2 



(128) 



diverges as (p2) « l-^ for Z2/A2 ^ 0. The limit I < A 
clearly demonstrates that (Pf) can detect properties of the 
probability distribution that are very difficult to be found 
by measuring powers of Xd and Pei- 

The expectation value {Pg) refiects properties of the 
probability distribution, but it cannot be expressed as a 
standard classical observable. It is not possible to assign 
some fixed value Pg{x,p) to every classical state {x,p), 
such that (p2^ — J^^w{x,p)Pg{x,p). Nevertheless, for a 
suitable class of probability distributions all moments (P") 
are calculable from w, such that Pg may be regarded as an 
observable in this sense. We will call observables of this 
type "statistical observables". They have a status similar 
to entropy or temperature, which can be assigned to equi- 
librium ensembles, but have no fixed value in a given mi- 
crostate. The question if p2 can be measured depends on 
the existence of measurement devices that can be brought 
into correlation with values of Pg in a given range. For- 
mally, we may associate to Pg the "nonlinear observable" 
{dxlnwy/i, with 



(P 



{dx InwY 



(129) 



x,p 



We note , however , that {Pf ) differs from 
l^^w{dx\nwy/16. 
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A very interesting feature of the use of the classical wave 
function 'ip{x,p) and the quantum formalism is the possibil- 
ity to express statistical observables as in the standard 
way as quantum operators. This allows us to use all the 
concepts from quantum mechanics, as the spectrum of the 
observables or the commutation relations and uncertainty 
relations. In particular, we can make a change of basis such 
that Ps and Xs are represented as diagonal operators. This 
defines new states of the system for which and Xg take 
fixed values. In the quantum formalism, Pg is a perfectly 
acceptable observable. If it can be measured, the predic- 
tions for such measurements are the ones from quantum 
formalism. 

We will see in sect. IVIIII that we can indeed define quan- 
tum operators for position and momentum 

Xq = Xci + 

Pq = PcI + Ips , [1q,Pq] (130) 

whose expectation values can be computed from the clas- 
sical wave function 'ip{x,p) and therefore from the classical 
probability w{x,p). In particular, the commutator obeys 
the standard quantum mechanical relation for the position 
and momentum operators. 

6. Classical wave function 

The relation w = ip'^ fixes iJj only up to a sign s{x,p) = 
±1, 

ip{x,p) = s{x,p)w^^^{x,p). (131) 

We will argue next that s{x,p) is severely resticted by 
properties of continuity and differentiability. Essentially, 
ip{x,p) is determined by w{x,p) and does not contain ad- 
ditional information. (Of course, the overall sign of ■)/' is 
arbitrary, but it plays no physical role.) If a continuous 
and arbitrarily often differentiable function tl){x,p) exists 
for a given choice of s{x, p) , it is clear that a different choice 
of s{x,p) would not share these properties any longer (ex- 
cept change of overall sign). 

We can discuss this issue in terms of the statistical ob- 
servables Xg, Ps or the "quantum observables" Xq, Pq. In- 
deed, the^ requirement that expectation values of the type 
{F{Xq, Pq)) can be computed in terms of tp{x,p) imposes 
restrictions on the sign s{x,p). If the expectation value 
(Pg) or (Pq) exists, this supposes that expressions of the 
type 

fidxi^)^ ^\ f w{dx\nw)^ (132) 

are well defined. Thus ip should be a continuous function 
of X, since a discontinuity would generate a divergent ex- 
pression on the l.h.s. of eq. (|132p . Similarly, the existence 
of 

Ifx^r ^ll^n^[dl\nw+^idx In w)^] ' (133) 

requires dxip to be continuous. Close to a zero of ip at xq 
this allows ip — a{x — xq), while a different sign for x < xq, 



as ip = a\x — xo\, is not consistent with eq. (|133p . We can 
define the expectation values of statistical observables as 
(Pg) , (Pg) in terms of w{x,p). If they are finite, an ap- 
propriate choice of s{x,p) exists such that the expressions 
of the expectation values in terms of operators acting on 
Tp^x^p) also exist. This restricts s{x,p) by the properties 
of w{x,p). 

If we assume that arbitrary powers (PJ') , {Xg) exist, the 
sign function s{x,p) is determined by the topology of the 
zeros of w{x,p). For any w{x,p) without zeros continuity 
of Ip requires that s{x,p) has to be the same for all points 
(x,p). Thus a sign flip of s{x,p) can only occur for zeros 
of w,w{xq,po) = 0. Assume next that w(x,p) has an "iso- 
lated zero" at {xo,Po) in the sense that w{x,p) is strictly 
positive in a neighborhood around (xo,po) (excluding the 
point {xo,po)). Then we may draw a circle around (xo,po), 
arbitrarily close to {xo,pq). On this circle w{x,p) is strictly 
positive, and continuity implies that s{x,p) cannot change 
sign on the circle. Extending its value to (xo,po) the sign of 
s{x, p) is the same in a whole neighborhood around (xp , Po) ■ 
This argument extends to zeros on a compact subspace of 
the phase space (with coordinates x,p). The sign must 
be the same in the region surrounding the subspace where 
w is strictly positive, and we can formally extend this sign 
to the subspace where w vanishes. 

Typically, w vanishes on the "boundary of phase space" 
for (a;^ + p'^) — >■ oo. We may have a line of zeros from 
one "point" (direction) on the boundary to another, divid- 
ing phase space into two pieces. At a given point (xq^pq) 
on the line we may denote by y a coordinate in the direc- 
tion perpendicular to the line, with |y| the distance from 
the line. The existence of (Pg) and (Xg) requires that 
^{dy/w / dy) vanishes for y — 0, while the existence of 
higher powers {Pf) etc. implies that the limits for y — > 0± 
of dy/w/dy exist. (The limits may be different for posi- 
tive and negative y.) For nonzero d^/w/dy the sign s{x,p) 
must jump at {xq,Po), such that ip has opposite sign on 
the two different sides of the line, limj,_>.o ^p — > ay. If the 
line is isolated (no crossing with other lines) it divides 
into two pieces with a jump of s at the division. This 
type of arguments can be extended to crossing lines and 
points on the line where a = \im.y^o{d'tp / dy) vanishes. For 
{d^/w/dy){y ^ 0) = and {d^ ^/w / dy"^ {y ^ 0) ^ there 
is no sign jump, with '\p ^ by^. The problem reduces to the 
sign of an infinitely often differentiable function 'ip{y) with 
a zero at y = 0. This type of discussion also extends to 
p-dimensional phase space, where the line is replaced by a 
(p — l)-dimensional hypersurface. 

In summary, for a large class of classical probability dis- 
tributions w{x,p) the classical wave function ip{x,p) can 
be determined in terms of w{x,p). The sign s{x,p) in eq. 
P3ip is computable from w{x,p). The expectation val- 
ues of non-commuting quantum operators, as Xq,Pq in 
eq. (|130l) . or more generally F{Xq, Pq), can be computed 
from the classical probability distribution. When expressed 
in terms of the classical wave function, they obey the stan- 
dard laws of quantum mechanics. 
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7. Interference for classical particles 

The particle- wave duality is not a characteristic of quan- 
tum particles. Any probabilistic description of particles 
induces the concept of a particle-wave. These wave as- 
pects are shared by classical particles as well. Even though 
particles can be viewed as discrete objects, the probabil- 
ity distribution w{x,p) is a continuous function or a field. 
The time evolution of w specifies the field equation for the 
wave. We also have introduced the wave function t{j{x,p) 
for classical particles. It has the properties of a probabil- 
ity amplitude, with w — The time evolution of w is 
mapped to a time evolution of "0 - it has the form of a wave 
equation or Schrodinger equation. 

If classical particles can be described by waves, with the 
same interpretation as probability amplitudes as in quan- 
tum mechanics, one may wonder what happens with pos- 
sible interference effects that are usually considered as be- 
ing characteristic for quantum particles. The generalized 
Schrodinger equation (|88p is linear in ip such that the su- 
perposition principle holds: if i/'i and are two solutions 
for the Schrodinger equation, also ipi + 'ip2 is a solution. If 
the Hamiltonian is given by the time evolution 

for the probabilities w = il>'^ is given by the linear Liou- 
ville equation, such that the superposition principle holds 
for the probabilities as well. One may ask if in a double 
slit experiment the probabilities add (no interference) or 
the amplitudes add (interference). We will see that this 
depends on the precise initial conditions for the classical 
probability distribution. 

Consider two initial classical wave packets wi{x,p,t = 
0) , W2{x^p,t — 0), chosen such that for wi the parti- 
cle passes through slit 1, and for W2 it passes through 
slit 2. This is possible since we may have wave packets 
which correspond to sufficiently focused beams (with ap- 
propriate small Aa;,Ap), such that the time evolution of 
wi{x,p,t) leads for all t to a vanishing probability at the 
location Xs2 of slit 2, wi{xs2iP,i) ~ 0. The initial wave 
functions corresponding to wi and W2 are ipi^x^Pyt = 0) 
and ip2{x^p,t = 0). If the initial conditions are set by 
w{t = 0) = wi{t = 0) or w{t = 0) = W2{t = 0) the particle 
will indeed "pass through only one slit" and no interference 
pattern will be obtained on a screen behind the slits. 

This situation extends to an initial probability distribu- 
tion 

w{t = 0) awi{t = 0) + (1 - a)w2{t = 0). (134) 

Due to the superposition principle for w the probability for 
finding hits on the screen at the location which correspond 
to particles which have passed through slit 1 amounts to a, 
and to 1 — a for the location corresponding to slit 2. The 
initial wave function for this setting is 

ip(t = 0) = {awi{t = 0) + (1 - a)w2it = 0))^^^ (135) 

One may, however, prepare also a different initial prob- 
ability distribution, 

V'(t = 0) = aii/)i(t = 0) + a2V^2(< = 0), 
w{t = 0) = ajwiit = 0) + alw2{t = 0) 

+2aia2Mt = 0)Mt = 0), (136) 



with a2 determined by ai and the normalization of tjj or 
w. Now the superposition principle for ijj implies tp{t) = 
aiV'i(t) + a2i'2{t) or 

w{x,p,t) — alwi{x,p,t) + alw2{x,p,t) 

+2aia2ipi{x,p,t)ip2ix,p,t). (137) 

We recognize an interference term ^ aia2 which leads, in 
principle, to an interference pattern on the screen. The 
interference can be positive or negative. It becomes an im- 
portant effect if ai and a2 are of similar magnitude and V'l 
and ■02 have a substantial overlap for some positions on the 
screen. The details of the interference differ, however, from 
the quantum particle. They also depend on the question 
which location observable is appropriate for measuring the 
location of the particle on the screen, e.g. Xci or Xq. 

VII. QUANTUM AND CLASSICAL PARTICLES 
IN A POTENTIAL 

In this section we compare the propagation of a quantum 
particle and a classical particle in a given potential V{x). 
We use for both a probabilistic description with a distribu- 
tion of initial values for location and momentum at time 
to = 0. We also employ for both the quantum mechanical 
formalism, keeping in mind that the latter can ultimately 
be obtained from a classical statistical ensemble. Using the 
same formalism constitutes an appropriate framework for 
a discussion of the differences between quantum and clas- 
sical particles. We have seen already that the expectation 
values of position and momentum obey the same evolution 
equation ([?T|) . Possible differences must then be connected 
to a different evolution of the probability distribution, or 
differences in the allowed initial conditions. 

1. One particle phase space distribution 

As a starting point for the classical particle we consider 
the probability in phase space w{x,p), from which observ- 
ables involving arbitrary functions of position and momen- 
tum can be computed 

{F{x,p)) ^ j I F{x,p)w{x,p). (138) 

J X J p 

(Here we use in d dimensions the shorthands — J d'^x, 
/p = (27r)~'' / d'^p. We keep our discussion one dimen- 
sional, but many formulae can be extended to arbitrary d 
if suitable scalar products of vectors are used when appro- 
priate.) 

We will compare w{x,p) with the Wigner representation 
Pw{z,q) of the density matrix for a quantum particle, as 
given for the free particle by eq. (IMl) . From the density 
matrix p{x, y) in position space the Wigner representation 
obtains as 

p^(z, q)^ j d{x - y)e-'''^-~y^p{x, y), (139) 

where z = {x + y)/2 is kept fixed for the integration over 
the "relative coordinate" {x ~ y). We next recall in simple 
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terms that appropriate correlation functions for a quantum 
particle obey eq. (|138p if w{x,p) is replaced by p^{x,p) 

The expectation value of an observable which is a func- 
tion Fx{x) of the position obeys 



{F,{x)) = tr(pF,(X)) = / p{x,x)F,{x) 

J X 

F,{x)p^{x,p). (140) 



X J p 



Here we use the inverse Fourier transform in order to ex- 
press p{x^ y) in terms of the Wigner representation p{z, q) 



Jp{x-y) 



Pu 



x + y 



with 



p{x,x) = / Pw{x,p), 
Jp 



,P , (141) 



(142) 



and X is the usual position operator in quantum mechan- 
ics. Similarly, we find for a function Fp(p) of momentum 

{Fp{p)} = tT{pFp{P)) = \hnJ FJ~i^]p{x,y) 



Fp{p)Pw{x,p), 



(143) 



X 'J p 



with P the quantum mechanical momentum operator. For 
the last equation (|143p we assume the absence of boundary 
terms, as appropriate if pw{x,p) has support only in a local 
region of space. Due to linearity, eqs. (11401) . (|143l) extend 
to all observables of the form F{x,p) = Fx{x) + Fp{p), and 
we find the expression analogous to eq. (|138|) 



(F(x,p)) 



F{x,p)p^(x,p). 



(144) 



For observables involving both powers of X and P the 
order of the operators matters. One finds 



tT{XPp) = J J x\^p- -d^] Pw{x,p) 



xp+ Pw{x,p), (145) 



while 



tv{PXp)^ J J i^xp^ -j pUx,p). (146) 

In particular, the correlation function for subsequent mea- 
surements of X and p [2\ , 



{XP)rn = {{X,P}p) 



xppw{x,p) (147) 



is given by the same expression as the correlation function 
for a classical particle if replaces w in eq. (I138p . The 
expectation values of arbitrary sequences of operators X 
and P can be represented as 

{S[X,P])= f f S[Xq,Pq]pUx,p), (148) 

where 

and the order in the sequence is the same on both sides 
of eq. (|148p . As it should be, the commutation relation is 
transferred from the operators X and P to Xq and Pq, 



[Xq,Pq] = 



(150) 



The r.h.s. of eq. (I148P assumes that partial integrations 
can be performed, i.e. that boundary terms where deriva- 
tive operators dp or dx stand on the first place in the se- 
quence ("on the left") vanish. 

In particular, we may consider the totally symmetrized 
functions of operators Fs{X, P), where we assume that an 
expansion of F{x,p) in powers of x and p is possible. It 
is defined by associating to each factor the totally sym- 
metrized combination of all possible ordered sequences, as 



(X^P^)s = -{X^P^ + P^X^ + XP^X + PX^P 
6 



'XPXP + PXPX). 



(151) 



Employing methods similar to [20] one finds the simple 
expression 



{FsiX,P)) 



F{x,p)pw{x,p). 



(152) 



This generalizes eq. (|144p . Comparison with eq. (|138p 
demonstrates the close formal correspondence between 
classical and quantum particles. We may define the one- 
particle phase space distribution fi{x,p) as fi{x,p) = 
w{x,p) for classical particles and fi{x,p) = pw{x,p) for 
quantum particles. The function fi{x,p) generates all ex- 
pectation values of totally symmetrized products of observ- 
ables according to eq. (I138P or p52l) . 

2. Differences between quantum and classical 
particles 

For a given distribution /i(x,p) the only difference be- 
tween a quantum particle and a classical particle concerns 
the ordering of sequences of measurements of the observ- 
able X or P. Only the totally symmetrized orderings coin- 
cide. For the classical particle the ordering does not matter 
for the classical correlation functions since Xd and Pd com- 
mute. This does not hold for a quantum particle, where 
the measurement correlations depend on the ordering of 
sequences of measurements j^j- For the quantum particle 
one finds for the measurement correlation for sequences of 
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measurements in different orders (the first measurement 
corresponds to the observable on the right of the sequence) 



(PXPX) 



{XPXP)m - {X^P^)ra - {P^ X^) m 

itr{p([X,P])'} = -i. (153) 



It is possible, however, to define also for the classical par- 
ticle a different set of correlation functions. We show in 
appendix B such an alternative definition which shares all 
properties of the quantum correlations. 

A second important difference between a quantum par- 
ticle and a classical particle concerns the time evolution of 
the one-particle phase space distribution. For the classical 
particle this is given by the Liouville equation 

dMx,P) = \--d. + ^^^pl w(x,v). (154) 

TO OX \ 



In contrast, for a quantum particle one finds 

dtPw{x,p) = \ -—dx +iV [x- l-dp 
TO V 2 



-iV \^+ ^dpj) pw{x,p) (155) 
^dx + '2V{x)s\TLi]^dxdp] \ pw{x,p). 



For this purpose we assume the von-Neumann equation for 
the time evolution of the density matrix for the quantum 
particle 



with 



dtp = -i[HQ,p\ 



(156) 



(157) 



This yields, for z = (.t + y)/2, the relation 
dtp{x,y) = / dtp^,{z,p)e'P^--y^ (158) 

tV{x) + tV{y) - LqX p^{z,p)e'P^^-y\ 
m J 



from which we infer eq. ()155p by partial integration. 

For a formal Taylor expansion of V in eq. ()155p in idp 
one finds 



Pw{x,p) 

(159) 



The first term agrees with the classical particle (|154p . In 
particular, for a free particle or for a harmonic potential 
V ^ x^ there is no difference between the time evolution of 
fi{x,p) for a classical or a quantum particle. For a given 
initial value of fi (x, p) the time evolution of all totally sym- 
metrized correlation functions is then the same for classical 
and quantum particles. 



In summary, the time evolution of the one particle dis- 
tribution for classical and quantum particles has a similar 
structure. For the classical particle it is given by the Liou- 
ville equation 



dtw = —iHLW, 



(160) 



with given by eq. (|89|) . The quantum particle is char- 
acterized by a different evolution operator 



dtPi, 



'iHwPw, 



with Hw obeying 



Hw = + V ix + ^dp j - V (x - ^dp 



(161) 



(162) 



Finally, a third difference concerns the allowed phase 
space distributions fi{x,p). In both cases the normaliza- 
tion is the same fi{x,p) — 1. For a classical parti- 
cle u) is a probability density in phase space and therefore 
fi{x,p) > must hold everywhere in phase space. For 
a quantum particle pw{x,p) is not anymore a probability 
density. Only the diagonal elements of the density matrix 
must be positive, i.e. 

pix,x) = / Pn,{x,p) > , pip,p) = / Pn,ix,p) > 0. 



(163) 

This does not exclude pw{x,p) to be negative in certain 
regions in phase space. 

On the other hand, not every classical fi{x,p) corre- 
sponds to a quantum mechanical density matrix. While 
the positivity and normalization of the diagonal elements 
p{x,x) or p{p,p) (eq. (|163p ) is obeyed automatically, this 
is not sufficient to guarantee the positivity of the density 
matrix. For example, a quantum mechanical density ma- 
trix has to obey the constraint 



tip' < 1 , 



Pix,y)p{y,x) 



\p{x,y)\'<l. (164) 



In terms of the Wigner representation this condition reads 

V-/ \p{x,y)\^^ f \pUz,p)\^<1. (165) 

Thus an ensemble of free classical particles can mimic the 
one particle phase space distribution for quantum particles 
with pw = w only if the "purity constraint" 



w^{z,p) < 1 



(166) 



z,p 



is obeyed. In this case one may expect the quantum me- 
chanical uncertainty relation for the location and momen- 
tum observables to hold. For a Gaussian initial wave packet 
the purity constraint reads (w^ = l-iAl*) 



(AojAp) exp <j — — )■ exp 



2A^Ap 



< 1. 



(167) 



21 



This is precisely the uncertainty relation A^Ap > 1/2, 
which is not automatically obeyed by all classical w^x^p). 

Indeed, a quantum mechanical density matrix has to 
obey a positivity constraint, namely that all eigenvalues 
of the matrix y) must be positive or zero. In turn, this 
imposes a constraint on the classical probability densities 
which reproduce a phase space description for a quantum 
particle: the hermitean matrix 

w{x, y)= J e'P^'=-y'>w (168) 

should only have positive or zero eigenvalues. We recall 
the properties 

zZ'(x,x) — / w{z — x,p)>0, 
Jp 

trw = / ■w{x,x) = 1. (169) 

Provided the classical probability distribution 'w{z,p) is 
square integrable and therefore obeys the purity constraint 
(|166l) . also the condition trw^ < 1 holds. If we consider the 
continuous function w^x, y) as the limit of a sequence with 
a finite number of degrees of freedom Af, the purity con- 
straint (|166l) follows from the positivity of w^x^y) by a 
suitable unitary transformation, 'w{x,y) — p{x)Sx,y The 
conditions p{x) > 0,J2xPi^) — ^ imply J^xP'^i^) — ^ ^'^'^ 
therefore ti w'^ < 1. Inversely, the constraint (|166p and 
the properties (|169p are not suflacient in order to guaran- 
tee the positivity of p{x,y), the latter being the strongest 
condition for fi(z,p) to describe a quantum particle. 

We conclude that the sets of allowed phase space distri- 
butions for quantum and classical particles have an overlap, 
namely whenever fi(z,p) > and the Fourier transform 
fi{x,y), given by eq. (|14ip or (|168p . is a positive ma- 
trix. Certain states are allowed only for classical particles 
(/i(z,p) > 0,/i(x, y) has negative eigenvalues) or only for 
quantum particles {fi{z,p) becomes negative in certain re- 
gions of phase space, /i(a;, y) is a positive matrix). Certain 
functions /i, where neither fi{x, y) is a positive matrix nor 
fi{z,p) is positive semidefinite, can be realized neither by 
a quantum nor by a classical particle. It would be inter- 
esting to know if generalizations of the particle concept 
beyond classical and quantum particles could realize such 
functions. 

Despite the close formal similarity between quantum and 
classical particles we summarize that differences arise on 
three levels, (i) The different status of w and pw implies 
different restrictions on the allowed states, in particular the 
allowed initial values. While w is a diagonal density ma- 
trix, is the Fourier transform of a density matrix with 
respect to only the relative coordinate x — y and there- 
fore not a density matrix by itself, (ii) For unharmonic 
potentials the different Hamiltonians generating the time 
evolution of classical and quantum particles are responsi- 
ble for a different time evolution of the phase-space density 
fi{x,p). (iii) The non-commutativity of the operators for 
a quantum particle entails a dependence of correlations on 
the sequence of measurements. 



VIII. QUANTUM PARTICLES FROM 
CLASSICAL PROBABILITIES 

The conceptual unification of quantum and classical par- 
ticles suggests that both could be based on a common de- 
scription in terms of a classical probability distribution in 
phase space. This possibility can also be regarded in the 
light of the limiting processes discussed in sect. IIVI The 
quantum particle requires in a sense less information than 
the classical particle. In eqs. (|56l) . (|57l) the quantum ob- 
servables use only a subspace. In a separate paper we will 
show in detail that such a description of a quantum particle 
moving in an arbitrary potential is indeed possible. 

In this short section we only sketch how all the prop- 
erties of a quantum particle can be described in terms of 
classical probabilities in phase space. For this purpose we 
use the classical wave function ip{x,p), which contains the 
same information as the classical probability distribution 
w{x,p) = ip'^(x,p). As compared to the classical particle, 
the quantum particle needs two modifications: (i) the use of 
quantum observables instead of classical observables, and 
(ii) a different time evolution of the classical probability 
density in phase space. 

Quantum observables Xq and Pq are defined by the op- 
erators (|149p . They obey the usual commutation relations 
for the quantum operators for position and momentum. 
The expectation value for an arbitrary sequence of such 
observables can be computed from the classical probability 
distribution by the usual quantum rule 

{F{Xq,Pq))^[ ^{x,p)F{XQ,PQ)i;{x,p). (170) 

We recall that the quantum observables involve the sta- 
tistical observables Xg and P^. If the ensemble is charac- 
terized by the classical probability distribution for states 
associated to the phase-space points {x,p), the quantum 
position and momentum do therefore not correspond to 
classical observables with a fixed value in every state. (If 
we want to realize Xq , Pq as classical observables we have 
to implement them in an ensemble with a larger set of 
states, as discussed in sect. IIIIl ) 

For a quantum particle, the time evolution of the classi- 
cal wave function (and the associated classical probability 
distribution in phase space) is given by 

idtilj{x,p) ^ Hwipix,p), (171) 

where we have replaced the Hamiltonian Hl by a modified 
Hamiltonian Hw given by eq. (|162l) . Eq. (|17ip should be 
regarded as the new fundamental equation describing the 
dynamics of particles. In consequence, the classical proba- 
bility density 'w{x,p) obeys a new non-linear time evolution 
equation instead of the Liouville equation. We can express 
the Wigner representation of the quantum mechanical den- 
sity matrix in terms of the classical wave function by 

Pw{x,p) = (172) 




ip{x + -,p + s)^p{x + + s') cos(sV — sr'). 
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It is straightforward to verify that with eq. (|171|) the time 
evolution of obeys the usual time evolution for the den- 
sity matrix of a quantum particle, with Hamiltonian Hq 
given by eq. (11571) . In other words, the time evolution 
of pw as defined by eq. (|172[) obeys eq. (|155p . We leave 
the proof that pw{x,p)^ as defined by eq. (|172|) . obeys the 
necessary positivity conditions for the associated density 
matrix to a separate paper. We will also show that the 
one particle distribution function for the quantum parti- 
cle can be understood as a coarse graining of the classical 
probability distribution in phase space. 

In summary, the expectation values of the quantum ob- 
servables Xq and Pq and all their correlation functions 
obey all the relations for a quantum particle in a poten- 
tial, including their time evolution. Starting at some initial 
time io with a classical probability distribution which cor- 
responds to a given pw{x,p), all quantum laws for a quan- 
tum particle in a potential are obeyed for all times, includ- 
ing characteristic phenomena as interference and tunneling. 
These correlation functions are the only thing measurable 
in this system - demonstrating that quantum mechanics 
can be described in terms of a classical probability distri- 
bution in phase space. 



IX. CONCLUSIONS AND DISCUSSION 

We have realized the description of a quantum parti- 
cle in a setting of a classical statistical ensemble with in- 
finitely many degrees of freedom. While the classical statis- 
tical ensemble contains information about the particle and 
its environment, the typical quantum mechanical features 
emerge if we concentrate on the statistical description of 
the subsystem for the particle. For the subsystem we deal 
with "incomplete statistics" for which joint probabilities 
can no longer be used for the prediction of outcomes of 
measurements of two observables. This is the origin for 
the representation of such observables by non-commuting 
operators in the associated quantum formalism. We have 
derived the quantum mechanical operator representation 
of the relevant position and momentum operators for a 
quantum particle explicitly, starting from a simple "one- 
bit particle" , generalizing to "Q-bit particles" , and finally 
taking the continuum limit Q — )■ oo in order to obtain posi- 
tion and momentum operators with a continuous spectrum. 
Heisenberg's law for the commutator of position and mo- 
mentum operators and the associated quantum mechanical 
uncertainty principle follow. 

In the course of this construction we have seen that al- 
ternative selections of observables, which lead in the con- 
tinuum limit to commuting position and momentum op- 
erators, are also possible. The same system of states and 
observables (the same Hilbert space) can therefore describe 
both quantum and classical particles. Quantum and classi- 
cal particles are characterized, however, by different Hamil- 
tonians. The appropriate Hamiltonian for real particles in 
nature can be (and has been) tested experimentally. For 
example, the double split experiment shows the interfer- 
ence pattern characteristic for a quantum particle. In con- 



trast, classical particles allow for states with arbitrarily 
sharp values of the location and momentum observables. If 
the probability distribution is sufficiently sharp, only parti- 
cles with trajectories passing by one of the splits will pass, 
even if the separation of the two splits is of the same order 
as the quantum mechanical wave length of the particles. 
This has clearly been falsified by experiment. 

The description of both quantum and classical parti- 
cles by a common formalism gives rise to the interesting 
question why nature prefers quantum particles. Interesting 
hints may come from an investigation of stability proper- 
ties of the quantum and classical wave functions that we 
have not yet investigated. We also note that the quantum 
particle is in a sense "minimal" since it uses less commut- 
ing operators as the classical particle. Furthermore, for the 
quantum particle the energy and the momentum operators 
are directly related to translations in time and space, in 
contrast to the classical particle. 

As a further aspect of the conceptual unification of quan- 
tum and classical particles we have developed a quantum 
formalism for classical particles. It is based on the classi- 
cal wave function in phase space ^{x,p)^ which equals the 
square root of the probability density up to a sign. Within 
this formalism the basic equation describing the dynamics 
of particles is given by a type of Schrodinger equation for ijj. 
For classical particles this replaces Newton's equations for 
trajectories - the latter only emerge as a particular case for 
infinitely sharp probability distributions. The use of the 
classical wave function allows for a simple description of 
new types of statistical observables that measure properties 
as the "roughness" of the probability distribution in posi- 
tion and momentum space. The concept of particle- wave 
duality also applies to classical particles. In particular, we 
have discussed interference effects for classical particles. 

Finally, our formulation of the time evolution of a classi- 
cal probability distribution w{x,p) in terms of a quantum 
mechanical wave function 'ip{x^p) may be used as a formal 
tool for the study of classical probability distributions. For 
example, we can immediately give a simple proof that dis- 
tributions of non-interacting particles, moving in an arbi- 
trary potential V{x), cannot reach the thermal equilibrium 
distribution as time goes to infinity. The unitarity of the 
time evolution of -0 leads to obstructions. We can construct 
an infinity of classical observables A{x,p) which commute 
with H and are therefore conserved. Such observables obey 
the differential equation 

^d:,A{x,p) = ^dpA{x,p). (173) 
m ox 

A simple family of solutions of eq. (|173[) is 

Aix,p) = f{E) , E = Vix) + 1^, (174) 

with an arbitrary function f(E). Thus arbitrary functions 
of the energy are conserved [l8|. An infinity of conserved 
quantities, which have to keep their initial values, contra- 
dicts the assumption that for large t a probability distribu- 
tion of thermal equilibrium is reached, where all expecta- 
tion values of observables can only depend on two param- 
eters, namely the temperature and the chemical potential. 
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We also find an infinity of static probability distributions 
in phase space. The corresponding static wave functions 
obey 



with 



p dV 



(175) 



Again, all functions of the form V'(-E') are static solutions. 

The observation that quantum mechanical features can 
obstruct the approach of classical probability distributions 
to thermal equilibrium distributions has been rnade earlier 
for classical field theories with interactions jlfl]. Similar 
obstructions for interacting particles exist and we believe 
that it is likely that obstructions beyond the conserved en- 
ergy distributions are present. This does not contradict the 
fact that many correlation functions approach the thermal 
limit for t — >■ oo, as demonstrated by numerical solutions in 
simple systems [l^. Of course, such obstructions can also 
be seen in the standard classical formalism. However, the 
quantum formalism makes them particularly transparent 
and perhaps more easy to approach. This could also hold 
for other features of classical probability distributions. 

Classical statistics and quantum mechanics are two sides 
of the same medal, rather than mutually exclusive con- 
cepts. This opens the potential for cross-fertilization be- 
tween the two formalisms, since they both describe the 
same physical reality. We do not know yet how far the 
practical use of this unification of concepts will reach. 



APPENDIX A: EXTENDED HAMILTONIAN FOR 
CLASSICAL PARTICLES 

The Hamiltonian Hl (jM]) commutes with the "classical 
Hamiltonian" 

2 

H,i^^ + V{x) , [Hl,H,i] - 0. (A.l) 
Zm 

The expectation value of Hd is therefore conserved and can 
be associated with the classical energy. As we have argued, 
Hci has no infiuence on the evolution of w{x,p). We can 
therefore define an extended Hamiltonian for the classical 
particle 



P 

2to 



V{x) 



i—d- 
m 



X 



iV'ix)dp. (A.2) 



The piece Hd does not infiuence the time evolution of 
the expectation value of any "diagonal observable" A{x,p), 
which remains purely dictated by H^. However, the use of 
H brings the quantum mechanical description of the clas- 
sical particle even closer to the usual setting. For example, 
a stationary state with fixed energy E obeys 

i/ci^ = Ei: , Hl^ = , H^Ij = E^. (A.3) 

The stationary wave function takes the form 

= u>"'"/'^(e) exp(— iet), 

2 

<x,p) = ^ + y{^)^ (A.4) 



ew^'^ie) = Ew^/^{e) , / / w{e) = 1 



(A.5) 



For a fixed energy w^/^(e) is a (5-type distribution with 
w^f^ie ^ E) = 0. We note that the wave function (fT2T|) 
remains an eigenvalue of Hl with eigenvalue zero also for 
more general w^^^{e). This implies H^ — H^ip for such 
states, and H plays the role of the classical energy. 

Even more generally, the Schrodinger equation dtip = 
—i{Hci + Hl)iP has always a solution 

ipix, p;t) = w^/'^{x,p;t)expli (^ + V{x)] t\ , (A.6) 



with 



2m 



(A.7) 



For this solution one has {Hl) = such that the average 
classical energy obeys 



'4}*H'4}. 



(A.8) 



Eq. (|A.6p is not the most general solution of the 
Schrodinger equation. Nevertheless, eq. (jA.7|) continues 
to hold even for the most general solution. This can be 
recovered formally for the extended Hamiltonian (jA.2p 

dtw = dt{r^)^rdt^ + [dtr)^ 

= -tr [{Hci + HL)i^] + iJ2)v/]^ 

= -i{r[Hci^] - [Hdr]^j) - i{r[HL^A + [Hl^]^) 

= -Lw, (A.9) 

where we use H*i = Hd , H*j^ = —Hl and the fact that iHl 
is a real first order differential operator, while ■ip*\Hd4'] — 
[HdV]^= [p^/(2m) + V{x)\'w. With 

dtw = 2w^/'^dtw^'^ = ~Lw = ~-2w^/'^Lw^''^ (A.IO) 



we verify eq. (|A.7I) . 

The most general solution of the Schrodinger equation 
can then be written in terms of ip{x,p; t) in eq. (jA.6p as 



(A.ll) 
(A.12) 



With 

9tV = dtipe^" -\- idtctij} 
= -i{Hd + HL)^ 
= -i[{Hd + HL)^y - i^[i?Le'"] 

we infer that a evolves according to 

dta — —iHLa — —La. (A. 13) 

According to eq. (I116P one finds 



(Hr 



w{La) 



wdta. 



(A.14) 
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Since the phase a is not observable by the use of "classical 
observables" A{x^p) we may use 5 = 0, which remains con- 
served by the time evolution. With the "initial condition" 
a(x,p; to) = the solution (|A.6|) . (|A.7|) holds for all t. 

We recall, however, that there is no need for an extension 
of the classical Hamiltonian according to eq. (|A.2p . In the 
main text we will continue to describe the time evolution 
for the classical particle by the Hamiltonian Hl. A real 
wave function ip = ±^/w, with a = 0,7r in eq. (jll4l) . 
is then a natural setting. In this case the energy is an 
independent observable E = H^i, not identical with the 
generator of the time evolution Hl. It commutes with 
and is therefore conserved. 



APPENDIX B: QUANTUM CORRELATIONS FOR 
CLASSICAL OBSERVABLES 

Let us consider a statistical ensemble for classical parti- 
cles for which the initial phase space distribution is chosen 
such that w{x,y) in eq. (|168l) is a positive matrix. For a 
free particle or a particle in a harmonic potential the phase 
space distribution /i(z,p) is then identical to the one of a 
quantum particle at all later times. The expectation values 
of all classical correlation functions are equal to the sym- 
metrized correlation functions for the quantum particle at 
all times. There seems to be no experimental way to dis- 
tinguish between a classical and a quantum particle in this 
case, suggesting that the classical and quantum particle are 
actually identical. 

Nevertheless, due to the non-commutativity of the posi- 
tion and momentum operators in quantum mechanics, one 
can also define correlation functions for an order of opera- 
tors differing from the totally symmetrized ordering (I15ip . 
While certain sequences of anticommutators yield the sym- 
metrized product, as 

^{{{P^X},P},X] = \{{P',X},X} = {X'P%, 

(B.l) 

(and similarly for X and P exchanged), other orderings 
differ, as 

i{p\X'}^iX^P%^^. (B.2) 

For a quantum system, the choice of the appropriate cor- 
relation function depends on the specific question. There 
is no unique definition of a correlation function for two 
powers of X and P. 

The same holds for a classical particle as well. For 
the special case where the classical probability distribu- 
tion equals the Wigner transform of an appropriate quan- 
tum mechanical density matrix, w{x,p) = p^{x^p), it 
is straightforward to construct correlation functions for 
the classical observables which differ from the classical or 
pointwise correlation function, in complete analogy to the 
quantum particle. In order to see this, we assume that 
w{x,y) in eq. (I168p is a positive matrix and formulate 
for this case the quantum mechanical concepts and laws 



directly in terms of the classical probability distribution 
w{z,p) or its partial Fourier transform 'w{x,y). In partic- 
ular, a pure state matrix ^ w is realized by a clas- 
sical probability distribution with the "folding property" 
= w or 

w{x,v)w{v,y) = w{x,y). (B.3) 

J V 

It can be described in the usual way by a complex quantum 
mechanical wave function ipq{x), with 

w{x,y) = ipQ{x)ip*Q{y). (B.4) 

(We use there the index Q in order to distinguish <pq(x) 
from the "classical wave function" ijjc{z,p) which obeys 

IV'C(2,P)P = W{z,p).) 

We recall that we can associate to an arbitrary classical 
observable A{z,p) a quantum mechanical operator A 

A{z,p)^ A = A,{X,P), (B.5) 

where A^ denotes the totally symmetrized product of the 
operators in a Taylor expansion (cf. eq. (|151I) V The oper- 
ators X,P act on w{x,y) or '^q{x) as 

Xw{x,y) = xw{x,y) , Pw{x,y) = -idxw{x,y), 
XipQ^x) = xlpq{x) , P(Pq{x) = -id^ifigix). (B.6) 

We have already shown the quantum laws for expectation 
values 

{A{z,p)) = / Aiz,p)wiz,p)=tT{Aw) (B.7) 

-/ z,p 

= / Aw{x,y)iy=^ = / ip*Q{x)AipQ{x), 

J X J X 

where the last relation holds for pure states. In contrast to 
the classical observables A{z,p) the corresponding quan- 
tum operators A do not commute. We can use them in or- 
der to derive Heisenberg's uncertainty relation for all clas- 
sical probability distributions that can be described by a 
positive matrix ?Zi(x, y) such that they obey the purity con- 
straint (jl66|) . This demonstrates in a simple way how the 
quantum formalism with non-commuting operators arises 
from a classical ensemble describing a probability distribu- 
tion of classical particles. Correlation functions different 
from the classical correlation function can now be defined 
by different orderings of the operators X and P. 

We close this appendix by the remark that a judgment 
if a given classical probability distribution w{x,p) can cor- 
respond to a quantum state may be often rather involved, 
since the positivity condition for w{x,y) may be hard to 
verify. The opposite way is much easier. In order to 
compute the classical probability distribution which corre- 
sponds to a given quantum state, it is sufficient to compute 
the Wigner representation of the density matrix puj{z,p) 
and to check if this is positive. 
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